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ABSTRACT 

We search 0.02 deg^ of the invariable plane for trans-Neptunian objects (TNOs) 25 AU or more distant 
using the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope. With 22 ksec 
per pointing, the search is > 50% complete for meoew < 29.2. Three new objects are discovered, 
the faintest with mean magnitude m — 28.3 (diameter « 25 km), which is 3 mag fainter than any 
previously well- measured Solar System body. Each new discovery is verified with a foUowup 18 ksec 
observation with the ACS, and the detection efficiency is verified with implanted objects. The three 
detections are a factor ~ 25 fewer than would be expected under extrapolation of the power-law 
differential sky density for brighter objects, E(m) = dN/drndfl cx 10"™, a « 0.63. Analysis of the 
ACS data and recent TNO surveys from the literature reveals departures from this power law at both 
the bright and faint ends. Division of the TNO sample by distance and inclination into "classical 
Kuiper belt" (CKB) and "Excited" samples reveals that E(to) differs for the two populations at 96% 
confidence, and both samples show departures from power-law behavior. A double power law I](m) 
adequately fits all data. Implications of these departures include the following. (1) The total mass of 
the "classical" Kuiper belt is w O.OIOM^, only a few times Pluto's mass, and is predominately in the 
form of ~ 100 km bodies (barring a secondary peak in the mass distribution at < 10 km sizes). The 
mass of Excited objects is perhaps a few times larger. (2) The Excited class has a shallower bright-end 
magnitude (and presumably size) distribution; the largest objects, including Pluto, comprise tens of 
percent of the total mass whereas the largest CKBOs are only ~ 2% of its mass. (3) The derived 
size distributions predict that the largest Excited body should be roughly the mass of Pluto, and the 
largest CKB body should have to^ w 20 — hence Pluto is feasibly considered to have originated from 
the same physical process as the Excited TNOs. (4) The observed deficit of small TNOs occurs in the 
size regime where present-day collisions are expected to be disruptive, suggesting extensive depletion 
by collisions. The Excited and CKB size distributions are qualitatively similar to some numerical 
models of growth and erosion, with both accretion and erosion appearing to have proceeded to more 
advanced stages in the Excited class than the CKB. (5) The lack of detections of distant TNOs implies 
that, if a mass of TNOs comparable to the CKB is present near the invariable plane beyond 50 AU, 
that distant population must be composed primarily of bodies smaller than « 40 km. (6) There are too 
few small CKBOs for this population to be the reservoir of Jupiter-family-comet precursors without 
a significant upturn in the population at diameters < 20 km. With optimistic model parameters and 
extrapolations, the Excited population could be the source reservoir. Implications of these discoveries 
for the formation and evolution of the outer Solar System are discussed. 
Subject headings: Kuiper Belt — solar system: formation 
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1. MOTIVATION 

The nebular hypothesis for th e formation of planetary 
systems is nearly 250 years old l)Kantill755|) and yet ob- 
servational support for the model is relatively recent. In 
the standard scenario, solids in the disk surrounding the 
protostar begin to coagulate into macroscopic objects, 
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which accrete to kilometer sizes. When the planetesimals 
become massive enough for gravitational focusing, run- 
away accretion begins. In the oligarchic growth phase, 
accretion is limited by excitations in the population in- 
duced by the largest few objects. In a protoplanetary 
disk, these largest planetesimals can reach a few M^, 
sufficient to trap the nebular gas, and rapid growth of 
gas giant(s) can ensue. The nebular gas is cleared by the 
stellar wind, and the remaining planetesimals are scat- 
tered away by the giant planets. 

Today we have many observat ions of dust and 
gas disks around young stars llO'Dell fc BeckwithI 
Il997t iBeckwith. Henning. & Nakagawal 1200(1) . evidence 
that supports the nebular hypothesis. Addition- 
ally, observations of dust disks around somewhat 
older stars suggest the presence of a population 
of du st-producing pla n etesim a ls in those s ystems , 
e.g. ISmith fc Terrild itl^^: IGreaves fit oiJ dTool : 
IKoerner. Sargent, fc Ostrofl lj2001() . Some of these dust 
disks exhibit structures that can perhaps be ascribed to 
embedded planetary systems ( Kuchncr & Holman 2003) . 
There is also now abundant evidence for the final stage 
of accretion — planet formation — as extrasolar giant plan- 
ets have bee n detected by radial velocity an d transit 
observations ijMarcv. Cochran. &: Mavoill2000j) . Though 
the basic idea of the nebular hypothesis remains intact, 
each new round of observations has led to fundamental 
changes in our view of planet formation. The presence 
of gas giants at < 1 AU, for example, was not well antic- 
ipated by theory, and migration is now recognized as an 
important process. 

It is unfortunate that direct observation of < 1000 km 
planetesimals in extrasolar systems is currently infeasible 
and likely to remain so for many decades. Such observa- 
tion would likely reveal further failures of imagination in 
our modeling of the planctcsimal phase. Fortunately, a 
portion of the Sun's planetesimal population is preserved 
for our examination in the region beyond Neptune, where 
growth time scales are longer, the accretion process ap- 
parently did not proceed to formation of planets, and 
the influence of the giant planets was not sufficient to 
remove all small bodies. Study of trans-Ncptunian ob- 
jects (TNOs) provides "ground truth" for models of the 
accretion, collisional erosion, and dynamical evolution of 
planetesimal populations. True to form, the TNO popu- 
lation only vaguely resembles the preconception of a dy- 
namically pristine planetesimal disk. With > 800 TNOs 
discovered between 1992 and the present, it is clear that 
the TNO population has several distinct dynamical com- 
ponents, all of which appear to have eccentricity and in- 
clination distributions that are too broad to be the undis- 
turbed remnants of the primordial population. The TNO 
population contains unmistakable signatures of interac- 
tions with Neptune and perhaps other massive bodies. 
With further study we can hope to understand the dy- 
namical history of this region. 

The physical properties of the TNOs, particularly the 
size distribution, are indicative of the accretion process. 
Observations to date are consistent with a distribution 
of diameters D that is a power law. dN/dD cx D'' , with 
g = 4.0 ± 0.5 l|Truiillo et all l200ll TJL). This distri- 
bution must fail at some I? > to avoid a divergence 
in the mass or reflected surface brightness of the trans- 
Neptunian cloud, but the scale of the breakdown in the 



power law is not usefu lly bounded by these constraints 
ijKenvon & WindhorstI 120011. In the current dynami- 
cal environment, TNO collisions are erosive for objects 
with diameters < 100 km, so that small objects have 
been removed from the population since the events or 
processes that excited the TNO dynamics (|Sterni.l99 5'). 
Rather soon after the discovery of the Kuiper Belt, there 
was speculation that th e size distribu tion might break at 
^ 50 km sizes ( Weissm an fc Levisonl [T997'l. but observa- 
tions to date have not evidenced this phenomenon. A 
generic prediction of accretion/erosion models is a break 
to a shallower size distribution below some size, but the 
size break is dependent upon fact ors such as the du ration 
of the accretion epoch (FarincUa . Davis, fc Sternl l2GG0'l. 
The mass in the trans-Neptunian region must have been 
substantially larger in the past in order to support the 
migration of Neptune (Hahn fc Malhotra 1999) and ac- 
cretion of the present TNO population (jStera il995() , 
but the relative importance of scattering and collisional 
grinding in mass removal is unknown. 

Extending our knowledge to the faintest (and hence 
smallest) possible TNOs is clearly desirable as there may 
be signatures of the collisional evolution or processes 
unanticipated by present theory. It is of further inter- 
est to see if the size distribution has dependence upon 
dynamical properties, as this can provide further insight 
into the dependence of the accretion/collision process 
upon the dynamics of the parent population. 

The Hubble Space Telescope (HST) is currently the 
observatory of choice for detection of the faintest pos- 
sible point sources. A detection of a very high density 
of mv > 27.8 TNOs using th e WFPC2 camera on HST 
is reported bv lCochra,n ^Hd] ((l995l) [CLSD] . For various 
reasons it is likely that these detection s were merely noise 
(jBrown. Kulkarni. & Liggett 1997; Glad man etaUil99^ 
(sec Bernstein. Allen, fc Trilling (2004) for further anal- 
ysis of the WFPC2 results). The installation of the Ad- 
vanced Camera for Surveys (ACS) on HST improved the 
field of view, efficiency, and sampling substantially. This 
paper describes the results of a large investment of HST 
time (125 orbits) into a search for TNOs using the ACS 
camera. 

Detection of faint objects requires long integration 
times, but a typical TNO moves the width of the HST 
point spread function (PSF) in only a few minutes. The 
ACS survey therefore uses a technique we call "digital 
tracking," in which a long series of exposures is acquired, 
with each individual exposure short enough to avoid 
trailing losses. The short exposures are shifted to follow a 
candidate TNO orbit and then summed, yielding an im- 
age with long exposure time that will detect TNOs on the 
chosen orbit with no trailing. The summation must be re- 
peated for all plausible TNO orbits that diverge by more 
than the PSF over the time span of the observations. 
This computationally intensive technique has been used 
successfully for seve r al ground-based faint- TNO searches 
(iTvson et all 11 992: Allen. Bernstein fc Malhotral 120011: 
iCladman et oiJil998: Chian g & Brown. 1999) and a vari- 
ant used by CLSD with WFPC2 data. We are able 
to detect TNOs to the fundamental limits set by pho- 
ton noise in the 22 ksec total exposure time of each 
ACS search field. The survey is > 50% complete for 
'71606 w < 29.2 mag, which is 2 mag fainter than any suc- 
cessful published TNO survey and 1.5 mag deeper than 
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the onset of false positives in the CLSD data. The area 
covered by the search is 0.02 deg^, 13 times the area of 
the CLSD search. The lessons learned from the ground- 
based and CLSD digital-tracking surveys have helped us 
to produce results that we believe are optimal and reli- 
able. 

The concepts and fundamental limits of digital track- 
ing in this and other applications are detailed in 
IBernstein. Allen, fc Trifl ing ( 2004). This paper summa- 
rizes the methodology of the ACS search, presents the de- 
tections and efficiencies, derives bounds on the apparent 
magnitude distribution of the TNOs and some dynamical 
subsamples, and discusses the implications for the evo- 
lution of the TNO system. Trilling & Bernstein ( 2004) 
present the variability data for the objects detected in 
the ACS survey. Ber nstein, Zach ar ias . fc Sm ith (2004} 
examine the current state of the art in astrometry for 
moving objects, and the utility of high-precision astrom- 
etry for orbit determination. 

2. DETECTION TECHNIQUES 

The search for moving objects to the photon-noise 
limit of a 22,000 s ACS integration requires a sophisti- 
cated analysis, attention to detail, approximately 30,000 
lines of code, and several CPU-years' worth of com- 
putation on 2.4 GHz Pentium processors. The unique 
tools of this data reduction are d escribed in detail in 
IBernstein. Allen, fc Trinind lj2004() . but we summarize 
here the aspects that are important for understanding 
the results. 

2.1. Observations 

The survey covers 6 slightly overlapping fields of view 
of the ACS camera. The spacecraft is oriented so that 
detector rows and columns are aligned to the local eclip- 
tic cardinal directions. The 6 pointings are arranged in 
a 2 X 3 mosaic, with the long axis in the ecliptic N-S di- 
rection. The southern two pointings are labeled A & B, 
the central two C & D, and the northern two E & F. The 
ACS pixel scale is nominally O'.'OSO, and nominal cover- 
age of the full mosaic FOV is « 400" x 600" = 0.019 deg^ 
The exposures at a given pointing are dithered by non- 
integral pixel steps, up to a few pixels, in order to im- 
prove the sampling of the static sky objects. The imaged 
field is not contiguous because our dithers do not span 
the gap between the two ACS CCDs. 

The field location was chosen subject to a number 
of criteria. The field center, 14h07m53?3 -11°21'38" 
(J2000), is only 3' from the invariable plane. The field 
trails Neptune by 99°, within the libration region for 
perihelia of TNOs in 2:1 and 3:2 resonance with Nep- 
tune (,,Malhotra„1996t iChiang fc Jordanll2002(l . A known 
TNO, 2000 FV53, is within pointing A for the full ob- 
serving period, allowing us to verify our navigation and 
orbital calculations. The field is placed and the obser- 
vations timed to minimize the loss of observing time to 
moonlight and South Atlantic Anomaly crossings, and 
to place the field 88° from opposition at the start of the 
observing sequence (see below). 

All exposures were taken through the F606W filter of 
the ACS using the Wide Field Camera (WFC). In the 
period UT 2003 January 26.014-31.341, which we call 
the "discovery epoch," 55 x 400 s exposures were taken at 



each of the six pointings.^ During 2003 February 05.835- 
09.703, the "recovery epoch," an additional 40 x 400 s 
exposures were taken at each pointing. The two sets 
of observations, 88-83° and 77-73° from opposition, are 
chosen to straddle the transition from prograde to ret- 
rograde motion for most TNOs. Hence any discovered 
objects have maximal chance of remaining in the mosaic 
FOV for the full 15-day duration of the HST observa- 
tions, and the image trailing due to apparent motion is 
minimized. 

Individual exposures are 340-410 s long, averaging 
400 s. Five exposures fit into a typical HST orbit, with 
fewer during radiation-impacted orbits. A set of 10 or 
15 exposures is taken during each HST visit to a given 
pointing. The pointings are visited in the pattern ABAB- 
CDCD-EFEF-ABAB-CDCD-EFEF during each of the 
two epochs of observation. So pointing C, for example, is 
sampled sporadically, at intervals as close as 8 minutes, 
over a time span of approximately 24 hours, during the 
first CDCD set of visits. Approximately 2 days later, the 
CDCD set of visits are repeated. Then w 7 days later, 
the cycle repeats for the recovery epoch. 

A few shorter exposures of the six pointings and of the 
outskirts of 47 Tucanae are taken in order to map the 
WFC point spread function (PSF) and provide astro- 
metric tie-ins. The performance of HST and ACS during 
the observations was nearly flawless. Comparison of the 
47 Tucanae images before and after the TNO observing 
cycle showed negligible change in the PSF, so we use a 
time-invariant (but spatially dependent) PSF map. 

2.2. Preprocessing and "Bright" Object Detection 

Once the data are placed in the HST archive, they are 
prepared for the moving-object search as follows: 

1. Bias removal, flat-fielding, and bad-pixel flagging 
are done by the HST "on-the-fly" processing. En- 
gineering keywords are checked for guiding errors 
or other problems. The uncertainty images are cor- 
rected for some errors in the STScI pipeline, and 
we create a weight image with the value w at each 
pixel being 1/cr^ (where a is the pixel's flux un- 
certainty). The weight is zeroed for defective and 
saturated pixels, and the data and weight images 
are changed into flux units. 

2. Objects in in dividual exposures are ca taloged using 
SExtractor ijBertin fc Arnoutslll99^ . 

3. The 47 Tucanae exposures are used to pro- 
duce a map of t he PSF for the WFC 
ijBernstein. Allen, fc Trifling 2004^1 . 

4. WFC distortion maps from STScI or from 
lAndersoiil 1)20031) are used to transform pixel po- 
sitions into a local tangent plane for each expo- 
sure, to accuracy « 10 mas; a translation plus lin- 
ear transformation are derived for each exposure 
to register all the cataloged objects onto a global 
tangent-plane coordinate system centered on the 
mosaic center. 

^ Exposure times varied slightly due to spacecraft constraints. 
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5. All exposures from the discovery epoch are com- 
bined into a deep image of the fixed sky. This tem- 
plate image has 0'.'025 pixels that are square (no 
distortions) on the global tangent plane, so that 
the PSF is now sampled near the Nyquist density. 
Because there are 55 exposures per pointing in the 
discovery epoch, each template pixel has 10 or more 
contributing images, and the template noise level 
is well below the individual exposures'. Sigma- 
clipping eliminates cosmic rays and bright moving 
objects from the template images. 

6. We interpolate the template image to the location 
of each pixel of each individual exposure. The in- 
terpolated template is then subtracted from each 
exposure. At pixels with very high flux (centers of 
bright stars and galaxies), we zero the weight image 
because the residuals to the template subtraction 
will rise above the noise. Note that the individual 
exposures have not been resampled in producing 
these "subtracted images." 

7. Artificial TNOs are added to the subtracted im- 
ages. One of us (MH) produces a list of objects 
with orbital elements and light-curve parameters 
selected at random from a chosen range. The po- 
sitions, magnitudes, and motions of these objects 
are calculated for each exposure. The position- 
dependent PSF is trailed for the motion and each 
artificial object added into the subtracted images, 
with appropriate Poisson noise in each pixel. One- 
third of the objects on the list are later revealed 
to the searchers (GMB and DET) for use in tuning 
the search algorithms. The searchers remain blind 
to the other two-thirds of the artificial objects until 
after a final TNO candidate list is produced. 

8. The subtracted images are searched for potential 
bright TNOs as follows. A PSF-matched, compen- 
sated filter is scanned across each subtracted image. 
Using the weight image, we can calculate the sig- 
nificance V (i.e., the signal-to-noise ratio) of each 
candidate point source peak in the subtracted im- 
age. All peaks with \v\ > 3.5 are noted and the 

of a fit to the PSF is calculated. Those that 
sufficiently resemble the PSF are recorded to a file 
of bright TNO candidates, to be examined later. 
Note that real TNOs will not fit the PSF precisely 
because of trailing, so our criterion for matching 
the PSF is kept loose, and the vast majority of 
candidates are cosmic rays. 

9. The subtracted images are "cleaned" in prepara- 
tion for the faint-object search as follows. Every 
pixel in the subtracted image that deviates by more 
than 5(7 from the mean sky level is fiagged. All 
weights are set to zero within a 2-pixel radius of 
each flagged pixel. This effectively masks all cos- 
mic rays and non-Gaussian noise in the subtracted 
images, which is extremely important for avoiding 
false-postive detections in the faint-object search. 
This process also masks bright TNOs and asteroids; 
the former have already been detected, however, in 
the previous step. 



10. A "flux image" is now created for each exposure. 
The fiux image is created on a regular grid in the 
global tangent-plane coordinates. Each such grid 
point is mapped back to a pixel position on the 
masked subtracted image, and we record the best- 
fit PSF fiux and its uncertainty for a point source at 
that location. Hence the "flux image" is a map of 
the brightness of a potential point source at any lo- 
cation in that exposure, and a weight (uncertainty) 
image is propagated as well. These flux images are 
the raw material for the faint-object search. 

Any potential bright moving objects must now be 
found on lists produced in Step|S| because the masking 
in Step ini may preclude their later detection. "Bright" 
in this context means detectable at > 3.5(7 in a single 
400 s HST exposure, which in practice corresponds to 
m < 27.6. We use here and henceforth the HST F606W 
magnitude system unless otherwise noted. The filter 
passband is roughly the union of V and R passbands, 
and the AB zeropoint is similar to a V zeropoint. 

Over 900,000 flux peaks trigger the f > 3.5 threshold 
in the discovery epoch. To fish the real (and implanted) 
TNOs from this sea of cosmic rays, we first require that 
a flux peak repeat in the same sky location (to 0'.'2) on 
successive exposures in one orbit, leaving 7700 pairs of 
detections. We reject linked detections that occur on 
the same detector pixels to avoid CCD defects. We next 
require two pairs of detections to exist within the same 
visit and be within « 2" hr~^ of each other, leaving 1300 
candidate quadruples of detections. Next a preliminary 
orbit is fit to each quadruple, and the methods of ti2.3.2l 
used to check whether the subtracted images are consis- 
tent with a point source moving on the putative orbit. 
This reduces the candidate list to 49 objects, of which 
46 are then revealed to be on the artificial-object list. 
The detection efficiency of the bright search for artificial 
objects is found to be 100% for m < 27.6. 

The 3 remaining objects are real: one is 2000 FV53, 
the previously known object, which at m = 23.4 is 
blindingly bright here, appearing at v ~ 80 in each of 
the 55 discovery epoch exposures and 40 recovery epoch 
exposures. The second bright detection is a new ob- 
ject, now given the preliminary designation 2003 BG91, 
with time-averaged magnitude (m) = 26.95 ± 0.02. The 
third detection from the bright search, 2003 BF91, has 
(m) = 28.15±0.04, but is highly variable and rises above 
the i> = 3.5 single-exposure threshold several times. 

The bright-object search is executed independently on 
the recovery epoch observations, revealing the same three 
objects, which are thus undoubtedly real. 

2.3. Faint Object Search 

The search for moving objects that are below the 
single-exposure detection threshold is much more com- 
putationally intensive. We must sum the available expo- 
sures along any potential TNO path through the discov- 
ery epoch exposures, then ask whether the best-fit flux 
for this path is safely above the expected noise level. 

2.3.1. The Search Space 

The space of TNO orbits is 6-dimensional, with one 
possible parameterization being {a, /?, d, d, /3, d}, where 
a and /3 are the angular position relative to the center 
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of the mosaic at some reference time Tq, c? is the geocen- 
tric distance at To, and the dots denote the TNO's space 
velocity in the same basis (c/. iBernstein fc Khushalanil 

((2OOO)). The hne-of-sight motion d has neghgible observ- 
able effect over the course of the 15-day HST observation, 
so we may set it to zero in our searches. This means we 
have 5 dimensions of TNO orbit space to search. We 
search on a grid of points in this space. The grid spac- 
ing in a and (3 is the pixel scale P of the flux images 
discussed above. The grid spacing Av in the velocity 
space (Qf,/3) should be fine enough that tracking errors 
are held to less than 1 pixel: Av < P/ AT, where AT is 
the time span of the observations being combined. Fi- 
nally we must choose a grid in distance d. The primary 
effect of d upon the apparent motion of the TNO is from 
the reflex of the Earth's orbit around the Sun (and the 
HST's orbit around Earth). The reflex motions scale as 
1 /d, so we choose a grid that is uniform in 7 = 1/d. We 
also note that the non-linear components of the TNO 
apparent motion all depend solely on 7 — primarily the 
reflex of Earth's orbital acceleration, but also the New- 
tonian gravitational acceleration of the TNO itself. The 
spacing A7 must be fine enough that errors in these non- 
linear motion components are held to <^ P. 

The number of grid points that must be searched then 
scales roughly as P~^AT^. We conduct our faint-object 
search in two passes: first, with P — O'.'OSO, and then a 
finer pass with P = O'.'OSO. The first pass runs quickly 
enough to be completed between receipt of the HST data 
in mid-February and scheduled foUowup observations at 
the Keck & Magellan telescopes in late April (see ij2.4() . 
But the PSF of the WFC is only 0'.'05 across, so mis- 
tracking by « 0.5P at P = O'.'OS causes significant blur- 
ring of the PSF in digitally tracked images, degrading 
our magnitude limit by « 0.2 mag. Hence we later run 
the finer-grid search to reach the ultimate limit of the 
WFC data. 

The bounds of the search space are determined as fol- 
lows: 

• We search 25 AU < d < 00. Even objects at d « 
1000 AU would move several ACS pixels over the 
course of our visits. 

• The perihelion of the orbit is constrained to be > 
10 AU. This places a lower limit on the transverse 
motion at a given d. 

• The orbit is assumed to be bound. This places an 
upper limit on the transverse motion at a given d. 

• The inclination of the orbit is assumed to be i < 
45°. This bounds the vertical component of the 
apparent motion. Note we search only prograde 
orbits. 

2.3.2. Steps for the Faint Search 

The faint search proceeds after Step^Jabove as follows 
for each the coarse P = 0'.'05 (discovery and recovery 
epochs) and the fine P = O'.'OS (discovery only) search: 

1. The flux images produced for this P in the search 
are split into six sets of visits. Set 1 contains the 
first ABAB sequence, set 2 the first CDCD se- 
quence, etc. The digital-tracking sums will be ac- 
cumulated over a set's worth of images, with time 



span AT < 24 hours. Digital tracking over the full 
5-day time span of the discovery epoch would be 
computationally infeasible. 

2. For each set, the outermost loop is over the dis- 
tance grid. The next inner loop is over the veloci- 
ties d and /3. At a given distance and velocity, we 
calculate an orbital shift for each exposure relative 
to the first exposure. The inner loops consist of 
summing the individual flux images at each pixel, 
with integer-pixel shifts defined by the velocity and 
distance. 

In the fine search, there are IS distance grid points 
and a total of ~ 7 x 10^ velocity grid points in the 
{a, $, d} space. For each set there are two pointings 
spanning « 1 x 10^ pixels in the flux images, with 
25-30 exposures per pixel per set. In total, the 
fine search tests « lO^'' points in the TNO phase 
space, requiring « 10^^ pixel additions to do so. 
This takes several CPU- years for 2.4 GHz Pentium 
4 processors, but a cluster of 10 CPUs at Penn and 
8 at Arizona reduces the required real time. 

3. At each grid point of the TNO search, the point- 
source fluxes along the track of the putative TNO 
from all exposures in the set are summed, as 
weighted by their inverse uncertainties, to form a 
total best-fit fiux and uncertainty. If the signifi- 
cance V = f /af exceeds a threshold of 4.0, the grid 
point is saved. 

4. Above-threshold grid points that abut in phase 
space are aggregated, and the most significant is 
saved. The output of the fine search is a list of 
K, 1.5 X lO'* significance peaks in the TNO phase 
space. 

5. For each detected peak a "tuneup" program is run, 
which fits a model moving point source to the pixel 
values in postage stamps from all subtracted im- 
ages. The gridded peak is the starting point and 
a, /3, d, and (3 are allowed to vary. The significance 
of real (or implanted) objects typically rises after 
tuneup since the optimized orbit is a better fit than 
the nearest grid point, and the position/velocity es- 
timates become more accurate. Significance peaks 
that are noise tend, however, to to become less sig- 
nificant, to have poor x^, and/or to fail to converge. 
The tuneup step reduces the number oi v > 4.0 
peaks in the fine search to 6 x 10^. 

6. The tuned-up peak catalog from one set is now 
compared to all other sets of the epoch; any pair 
of peaks that might correspond to a common or- 
bit are linked and passed to the next step. Note 
that TNOs which cross the boundaries of the ACS 
pointings are found as efficiently as those which do 
not. There are S x 10^ (non-unique) linked peak 
pairs in the fine search, of which S x 10** have total 
significance v >1 . 

7. All of the linked pairs with v > 7 are again run 
through the tuneup program, but this time all of 
the exposures from the entire epoch are used. The 
arc is now sufficiently long (typically S days) that 
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we can allow the distance d to vary without fear 
of degeneracy. A few detection candidates with 
— DOF > 150 in the fit of the moving-source 
model to the data are rejected; inspection shows 
these to be spurious detections near the residuals 
of diffractions spikes of bright stars. We apply a 
threshold of > 8.2 (i/ > 10 for the coarse search) 
to obtain the KBO candidate list of « 100 objects. 

The histogram of detections vs significance rises very 
rapidly below ly = 8.2, which is to be expected from 
Gaussian noise in a search of 10^^ or so p hase space loca- 
tions (,Bcrnstcin , Allen, & TriUin3 l200^ . The threshold 
is placed at the tail of this false-positive distribution. De- 
tection candidates above this threshold are inspected by 
eye, with 2-3 being clearly associated with subtraction 
residuals and other data flaws. 

In the coarse search, there are 92 detections with 
ly > 10. The bhnd list of artificial TNOs is then re- 
vealed, and 89 of 92 detections are found to be implanted 
objects. Two of the remaining detections, 2000 FV53 and 
2003 BFgi, coincide with bright detections. The last is 
a new object, 2003 BHgi, discovered with significance 
ly = 16.7 and mean magnitude m — 28.35. 

For the fine search (which had an independent set of 
implanted objects), there are 67 detections, of which 64 
are found to be on the list of implanted KBOs. The three 
remaining detections are again 2000 FV53, 2003 BF91, 
and 2003 BH91. 

The faint search technique is also applied to the re- 
covery epoch with a coarse {P = 0'.'05) grid. The same 
candidates are independently detected above the = 10 
threshold. 

Figure shows postage-stamp images of 2000 FV53 
and the three new detections, as we improve the depth 
of images by summing more exposures. 

2.4. Orbit Determination and Recovery 

Each of the TNOs detected in the discovery epoch is 
also clearly detected in the recovery epoch. We now 
combine the information from all exposures in the en- 
tire ACS campaign to get the best possible constraint 
on each object's orbit. We again invoke the tuneup pro- 
gram, whereby the orbital parameters are varied to maxi- 
mize the significance of the detection of the moving point 
source. More specifically, the orbital parameters deter- 
mine the location of the PSF and the degree of trailing in 
each individual exposure. The two endpoints of the trail 
are converted into pixel coordinates using the registra- 
tion information and the distortion maps. We calculate 
the PSF at the TNO location using the spatially vary- 
ing PSF maps from 47 Tucanae, and we smear this PSF 
to the required trail length. The flux of the TNO is al- 
lowed to vary in a stepwise fashion from orbit to orbit (or 
from exposure to exposure for the high-S/N 2000 FV53) 
A model with constant fiux for a given TNO would be 
a poor fit, as all of the detected objects have significant 
fiux variations. The moving, variable-fiux model is then 
fit to the subtracted images, with all orbital elements and 
fiuxes being optimized. A byproduct of this orbital op- 
timization is an optimally measured light curve for each 
object. Analvs is and int erpretation of these light curves 
is in .Trifling fc Bernstein (,2004) . 

For the final orbit determination, all six orbital pa- 



rameters are allowed to vary. The 2000 FV53 data are 
of such high quality — positional accuracy of ~ 1 mas for 
each of the 95 exposures — that the line-of-sight veloc- 
ity, and hence a a nd e, are significantly constrained with 
only a 13-day arc. iBernstein. Zacharias. fc Smith' (1 20041 ) 
will consider in detail the techniques, limitations, and 
beneflts of such high-precision astrometry for the deter- 
mination of Solar System orbits. 

For the three newly detected objects, the line-of-sight 
motion is still poorly determined over the 13-day arc. In 
the final orbit fit to the HST data, we include a prior con- 
straint on the kinetic and potential energies that weakly 
pushes the orbit to circularity: 

X^rior = 4(2KE/PE+l)2. (1) 

An unbound or plunging orbit is thus penalized as a 2a 
deviation. The result of the fitting process are best- 
fitting orbital parameters (in the {a, (3,...} basis) for 
each object and covariance matrices for each, which can 
be used as described in Be rnstein fc Khushalani (2Q(1Q3 
to give orbital elements and position pre/postdictions 
with associated uncertainties. 

We attempted retrieval of all new objects using the 
imaging mode of the DEIMOS instrument on the Keck II 
telescope on the nights of 27 and 28 April 2003. The error 
ellipses for all three objects fit within a single DEIMOS 
field of view, so for each object we have 5 hours of inte- 
gration in the R band on each of 2 nights. We sum the 
Keck exposures to follow the motion vectors predicted 
for the TNOs by the HST data. 2003 BG91 is detected 
at i? « 27 on each night, but 2003 BF91 and 2003 BHgi 
remain below the detection threshold. 

Attempts to retrieve the objects with the Magellan II 
telescope on 1-2 June 2003 were foiled by poor weather. 

The orbital constraints are now refined using the Keck 
position. Table Ogives the discovery circumstances and 
best-fit orbital elements for each object. They all have or- 
bits consistent with "classical Kuiper belt" objects (CK- 
BOs), with distances of 40-43 AU and inclinations of 
< 3°. The orbital eccentricities either are (2003 BG91) 
or are consistent with (2003 BF91, 2003 BH91) e < 0.08. 
It is interesting to note that no Plutinos were discovered 
despite the fact that our observations were in the lon- 
gitude region where Plutinos reach perihelion. Likewise 
no high-eccentricity or distant objects were found. The 
implications of their absence are discussed below. 

2.5. Detection Efficiencies 

We next address the important issue of whether our 
TNO search is complete {i.e. free of false negatives) and 
reliable (free of false positives.) 

2.5.1. Reliability 

We are claiming to have examined > 10^** possible 
TNO sites in phase space and have exactly zero false 
positives with v > 8.2. This is not a trivial issue, as 
more than one publication claiming detection of TNOs 
at R > 26 has upon furth er exami nation been found 
iBrown. Kulkarni. Liggett^l 997: Gla.dman et »ijl2001^ 
to have primarily false positive detections. In the ACS 
program, however, the detections are unambiguous: each 
of the 3 new objects is independently detected in the re- 
covery epoch as well as the discovery epoch. Further- 
more, the one sufficiently bright object is recovered at 
Keck. 
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Fig. 1. — Postage-stamp digitally-tracked images of all four TNOs detected in the ACS data. Successive rows show images 
with more contributing integration time, starting with a randomly selected single exposure and ending with the summed image 
of all available ACS data. All images are shown with the same greyscale. Lowest row gives final signal-to-noise ratio of each 
object, save 2000 FV53, for which the full-survey sum is omitted and the S/N per exposure is listed. Note that the faintest 
object is undetectable on single exposures, and yet is 0.8 mag brighter than our estimated completeness limit. 



Table 1. Properties and Barycentric Elements of Objects Detected 



Name 




a 


e 


i 


Mean F606W mag 


Diameter*^ 


(AU) 


(AU) 






(degrees) 


(km) 


2000 FVss'^ 


32.92 ±0.00 


39.02 ±0.02 


0.156 ± 0.001 


17.35 ±0.00 


23.41 ±0.01 


166 


2003 BG91 


40.26 ± 0.00 


43.29 ±0.06 


0.071 ± 0.004 


2.46 ± 0.00 


26.95 ±0.02 


44 


2003 BF91 


42.14 ±0.01 


50 ±20 


0.4 ±0.4 


1.49 ±0.01 


28.15 ±0.04 


28 


2003 BH91 


42.55 ±0.02 


45± 13 


0.2 ± 0.7 


1.97 ±0.02 


28.38 ±0.05 


25 



''Heliocentric distance at discovery. 

'^Assuming spherical body with geometric albedo of 0.04. 

"^Previously known TNO targeted for this study. Elements reported here are from ACS data alone. 



2.5.2. Completeness 

Is the search complete? This issue is addressed pri- 
marily through the implantation and blind retrieval of 
artificial TNOs. We implant two distinct sets of arti- 
ficial TNOs into the discovery epoch data: one for the 
coarse search and one for the faint search. In each case, 
the artificial TNO orbital elements are chosen at ran- 
dom from a constrained range of the element space. The 
range of elements is carefully chosen so that the artificial 



objects overfill the ranges of position, velocity, distance, 
and magnitude to which the search is sensitive.'^ From 
the randomly selected elements, we can then calculate 
the geometric search area by noting which objects fall 
into the field of view for the requisite number of expo- 
sures. From the final object list in each search, we cal- 
culate the probability of detection for objects that meet 

^ The inclination range of artificial TNOs is limited due to a 
software bug. We also do not place artificial TNOs at d > 200 AU. 
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Fig. 2. — Probability of detection vs mean measured magni- 
tude for the artificial TNOs in the fine search of the discovery 
epoch. The histogram gives the results from the lOf artificial 
objects within the FOV, and the curve is Equation 



the geometric criteria. The product of these two is the 
effective area QcS, which wiU be a function of apparent 
magnitude, and could depend upon such quantities as 
distance, rate of motion, and light-curve amplitude. 

For the coarse search, w 150 implanted TNOs have 
27.3 < m < 29.4 and light-curve amplitudes up to 
0.2 mag. In the bright search, 46 of these are recovered, 
and 89 are recovered in the faint/coarse search. From 
these we verify that there is no gap between the mag- 
nitude ranges for which the bright search and the faint 
search are 100% effective. The area lost to bright stars 
and galaxies is negligible because the PSF of ACS images 
is very small, and is also stable, so the fixed-sky subtrac- 
tion is very successful. The effective area has no detected 
dependence upon TNO distance or velocity within our 
TNO phase space search grid. This is expected since all 
TNOs should move several pixels from orbit to orbit, yet 
have average trailing loss < 0.1 mag. 

For the faint search, artificial TNOs are generated with 
28.6 < m < 29.4 in order to more carefully probe the 
limiting magnitude of the survey. A TNO is considered 
to be in the survey area if it is imaged in at least 3 of the 
HST visits of the discovery epoch; 101 artificial TNOs 
meet this criterion, of which 64 are detected ai v > 8.2. 
Figure 12 plots the recovery efficiency vs mean magnitude 
for the faint search. The effective area vs magnitude is 
well described by 

rJeff = (r!o/2)erfc[(m - TO5o)/2w], (2) 

where f^o = 0.019 deg^ is the peak effective solid angle, 
TO50 = 29.17 is the F606W mag at which the effective 
area drops 50%, and w — 0.08 mag is a transition width. 
The detection efficiency again has no measurable depen- 
dence upon distance or velocity over the search space. 

Why should we trust the artificial-KBO tests to verify 
our completeness? After all, if the implantation process 
and the search/extraction process make common errors 
in flux scale, orbit calculation, PSF shape, or image dis- 



tortion, then the artificial objects could be detected at 
high efficiency while real objects are not. 

We note first that the orbit calculation code used for 
object implantation was written by one of us (MH) while 
that used for extraction was independently written by 
another (GMB), and both codes were checked against 
each other and the JPL online Horizons service.^ 

The targeted TNO 2000 FV53 helps us to address con- 
cerns about errors in orbit calculations, image registra- 
tion and distortion, or PSF estimation. The individual 
exposures for 2000 FV53 are fit by our moving-PSF model 
to good precision: its positions match the extrapolation 
of previous observations to the accuracy of the extrapo- 
lation, and we find the positions consistent with a refined 
orbit to the level of « 3 mas, or < 0.1 pixel on the WFC. 
This is better than the claimed accuracy of the distortion 
map. We are thus reassured that our models and code for 
spacecraft navigation, image registration, orbital motion, 
and field distortion are correct to the accuracy required 
for the search. The images of 2000 FV53 are formally 
inconsistent with the PSF model (x^ per DOF is > 1) 
but this is because the S/N ratio of the 2000 FV53 obser- 
vations are very high. The deviations from the fit are at 
the level of a few percent of the PSF, and hence the PSF 
model is sufficiently accurate for the fainter detections. 

2.5.3. A Caveat on Light Curves 

The selection function for TNOs with variable mag- 
nitude is complex: the object must be seen during at 
least 3, preferably 4, HST visits with a S/N ratio of > 4, 
to survive the detection cuts. The timing of these vis- 
its is irregular, hence there is no simple way to quantify 
the impact of light-curve variation on detectability. The 
implanted TNOs were given sinusoidal light curves with 
peak-to-peak amplitudes chosen uniformly between and 
0.2 mag and periods chosen uniformly between 0.05 and 
1.3 days — in this range, we did not note any change in 
detection probability vs magnitude. We know, however, 
that there exist TNOs with light c urve amplitudes near 
or above 1 mag, such as 2003 BF91 fTrilling fc Ber nsteinI 
I2004|) and 2001 QG298 (ShcDoard & Jcwitt 2004), and 
we should investigate the effects of high variability on 
the detection properties of our and other surveys. Subtle 
biases on light-curve shape are present for all TNO sur- 
veys, though other authors have chosen, like us, to ignore 
them for simplicity. In Appendix IXI we demonstrate that 
these biases are too small to be significant with current 
data, but may be important for future larger surveys. 

3. CONSTRAINTS ON THE TRANS-NEPTUNIAN 
POPULATION 

The ACS survey detects a total of three objects (not 
counting the targeted 2000 FV53), described in Tabled 
over an effective search area described by Equation (01. 
TJL fit a power law to the cumulative ecliptic surface- 
density distribution of TNOs: 

A^(< R) = \{)°'(^-^-o) deg-^ (3) 

with R being the i?-band apparent magnitude, a = 
0.63±0.06 and Rq = 23.0. Their fit is to survey data from 
19 < i? < 27. Taking our limit = 29.17 to be equiv- 
alent to a limit of i? < 28.8 ( tl3.1|l . an extrapolation pre- 
dicts « 85 detectable objects in our survey. This is quite 



|http: //ssd. jpl ■ nasa.gov/horizons .html 
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inconsistent with our observation — even the TJL 2cr hmit 
of a = 0.51 predicts « 16 detections in our survey — and 
it is immediately obvious that the magnitude (and hence 
size) distribution of TNOs changes behavior somewhere 
in the 25 < to < 29 range. In this section we quantify 
the nature of this breakdown in the single power law, 
and calculate the implications for integral properties of 
the TNO population. 

3.1. Compendium of TNO Survey Data 

We wish to derive the differential surface density 
^{R) = dN/dRdO. of TNOs per R magnitude interval 
using all possible reliable published survey data. The 
requirements for published survey data to be useful are: 

• The coordinates of all fields searched must be given. 

• The effective search area as a function of m for each 
field must be given, preferably derived by Monte- 
Carlo tests. 

• The circumstances of discovery of all detected ob- 
jects must be given, including apparent magnitude, 
estimated heliocentric distance d, and estimated in- 
clination i of the orbit. 

Note that we do not require that all detected objects have 
fully determined orbits. While the R < 24 TNO discov- 
eries have been recovered with admirable completeness, 
the practical difficulties of recovery for fainter objects 
have precluded observational arcs longer than w 1 day 
for any object with R > 25.6 (prior to this ACS survey). 
Fits to one-day arcs yield d and i to 10-20% accuracy but 
other orbital properties are highly degenerate. Hence a 
comprehensive study of both and bright and faint detec- 
tions can as yet make use only of d and i to categorize 
the TNOs. 

The sky-plane density of TNOs is certainly a func- 
tion of latitude relative to the midplane of the popu- 
lation. Most of the R > 22 searches have been tar- 
geted to the ecliptic plane, or less frequently the in- 
variable plane. Brighter surveys cover larger area and 
have ranged farther from the ecliptic. Proper compar- 
ison of bright to faint TNO densities requires that we 
consider sky densities measured within a fixed band of 
TNO latitude. If the latitude distribution of TNOs were 
well known we could make use of all the available survey 
data. While the midplane of the TNO population has 
recently been estimated to be » .7 ± 0.4° from the in- 
variable plane ijBrown fc Panll200^ . the full distribution 
remains poorly constrained. We restrict the published 
survey data compendium to fields with invariable lati- 
tude < 3°. The TJL estimate of the inclination distribu- 
tion of bright CKBOs implies a drop by factor « 2 in the 
sky plane density from zero to three degrees ecliptic lati- 
tude, so there may remain substantial inhomogeneity in 
comparing surveys over a ±3° swath. Attempts to select 
a narrower latitude range are counterproductive, how- 
ever, given our poor knowledge of the sky distribution of 
TNOs. 

The resonant TNO population has longitudinal struc- 
ture in the sky-plane density as well, with more ob- 
jects being found at bright magnitudes in the direc- 
tions perpendicular to Neptune. Plutinos (3:2 Neptune 
resonators) have perihelion positions that librate about 



these points. This longitude variation has yet to be 
mapped in any way, so a correction is not possible. The 
effect upon our results is not likely to be significant, be- 
cause surveys at all magnitudes span a range of longi- 
tudes. The exception is our uniquely deep ACS field, 
which though pointed in the region of Plutino perihelion 
libration does not detect any Plutinos. Our most precise 
analyses will in any case be done on samples intended to 
exclude Plutinos. 

A subtle difference between the bright and faint TNO 
samples is that the former are typically discovered on 
single short (< 10 minutes) exposures, and hence mea- 
sure the instantanenous magnitude distribution. Objects 
with i? > 25 are detected in summations of many hours' 
worth of exposures, and depend upon the ffux of the 
TNOs averaged over their light curves (see also discus- 
sion m 12331). In Appendix m we show that this effect 
is insignificant for the current data. 

We list in Table El the published TNO surveys that 
meet the requirements. We restrict our consideration to 
those works that dominate the surveyed area at a given 
magnitude, and we omit surveys that have been shown to 
contain significant false-positive contamination. For the 
purposes of the S(i?) analysis, we wish to standardize 
all magnitudes to the R band. The La and CB data are 
reported in V band; Tcgler & Romanishin ( 2003f) present 
accurate colors for many (bright-end) TNOs, and the 
mean V — R is ~ 0.6 mag, which we apply to the V 
detections. Some of the ABM fields use a "V7?" filter, so 
for these we apply the color correction given by ABM, 
assuming again V — R — 0.6 for an average TNO. The 
F606W filter on the ACS WFC is essentially the union 
of the V and R passbands. iTegler fc RomanishinI 1)20031 ^ 
show that the average TNO is « 0.39 mag redder in V—R 
than the Sun, so we presume that the average mgoe — R 
TNO color is about 0.20 mag redder than Solar. Taking 
the Solar V — mQQQ = 0.06 mag from the ACS Instrument 
Handbook, an average TNO should have meoe — -R ~ 
0.4 mag. We correct the ACS limits and detections to 
R using this value. Henceforth we will use only i?-band 
magnitudes. In Appendix^we show that variance in V— 
R colors of TNOs has negligible effect upon our analysis 
of the current data. 

Some other adjustments to the published survey data 
are necessary: 

• The e ffective search area of each iLarsen et a,l\ 
l)200lD field is taken to be the product of its geomet- 
ric area and the F{T) entry denoting the fraction 
of the field that is estimated to be unique to the 
survey in their Table 1. We crudely fit a complete- 
ness model of the form (O to the completeness for 
each seeing bin listed in their Table 3. The effective 
area of all search fields centered within ±3° invari- 
able latitude are summed to give a total useful sur- 
vey area for the survey, and we only count objects 
detected in these low-latitude fields. The redun- 
dancy and broad latitude coverage of this survey 
mean that its peak effective area, for our purposes, 
is only 20% of its raw angular coverage. 

• In the TJL data, no distance or inclination infor- 
mation is available for 7 of 74 objects detected near 
the ecliptic plane. TJL note that this information 
is missing because of inclement weather at foUowup 
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time and therefore these 7 objects should be drawn 
from the same distribution as the remainder. We 
therefore omit these 7 from our hsting and decrease 
the tabulated effective areas by 7/74 = 9.5% in or- 
der to reflect this followup inefficiency. 



Allen. Bgrng^gii^^ Malhotr 
Allen. Bernstein fc Malhotr; 



are merged for this analysis. 



and 
[ABM] 



' iGladman et all l|200lD describes two searches, one 
with CFHT and one with the VLT. We sum their 
effective areas and detections in this analysis. 



• iTruiillo fc Brownl l|2003|) do not give individual 
field coordinates, but do give the total sky coverage 
as a function of invariable latitude, which suffices 
for our purpose. This preliminary report does not 
include a detection efhciency analysis, merely an 
estimate of a 50% completeness level R « 20.7. 
We avoid this uncertainty by making use of the 
TB data only for R < 20.2, and assuming that 
in this range the detection efficiency is a constant 
85% over the surveyed area. Note that the effec- 
tive search area of TB comprises the majority of 
the ±3° latitude region. 



• The brighte st surveys l|Truiillo fc Brownl l200l 

ILarsen et a/.ll200'l|l have inaccurate magnitudes for 
their detections, due to varying observing condi- 
tions and ill-defined passbands. Nearly all of these 
objects have, however, been carefully reobserved 
by other authors for color and variability informa- 
tion, and we can replace the original survey mag- 
nitudes with highly accurate i?-band mean mag- 
nitudes. The original magnitude uncertainties re- 
main relevant, however, for treatment of incom- 
pleteness, as discussed in Appendix 1X1 

• We truncate all the efficiency functions rj to zero 
when they drop below « 15% of the peak value 
for that survey, and ignore detections faintward of 
this point. In this way we avoid making our like- 
lihoods sensitive to rare detections in the (poorly 
determined) tails of the detection fimction. 

We define three dynamical groupings of the detected 
TNOs in these surveys: 

1. The TNO sample holds all objects discovered at 
heliocentric distances d > 25 AU. One known Cen- 
taur (1995 SN55) sneaks into this TNO sample, but 
we do not omit it because similar objects found in 
the faint sample would not have been rejected. 

2. The CKBO sample is the subset of the TNO sam- 
ple having 38 < d < 55 AU and i < 5° . This is 
intended to exclude resonant and scattered objects 
to the extent possible with our limited orbital in- 
formation. 

3. The Excited sample is the complement of the 
CKBO sample in the TNO sample. High inclina- 
tions and/or proximity to Neptune would indicate 
substantial past interactions with Neptune or an- 
other massive body. 



Note that we have used ecliptic inclinations rather than 
invariable, since the latter are not generally available. 
We have used the central values for d and i even when 
the surveys report uncertainty ranges that cross our def- 
initional boundaries. We have ignored the possibilities 
of overlaps in survey areas, and omitted targeted objects 
such as 2000 FV53. 

The TNO sample under analysis thus contains 129 de- 
tections spanning 19.5 < R < 28.0, of which 69 are as- 
signed to the CKBO class. Figure |3| shows the ilcff of the 
published surveys vs magnitude, and the binned magni- 
tude distribution of the detections. 

Our definition of the CKBO class is imperfect because 
we are restricted to use of d and i in classification. Res- 
onant and "scattered-disk" TNOs can also slip into the 
CKBO category under some conditions, and Centaurs 
near aphelion may be accepted as either CKBOs or Ex- 
cited TNOs. Of the objects classed as CKBOs in this 
study, 39 have sufficiently long arcs to determine a and 
e. Of these, all have 42 < a < 48 and e < 0.2 except the 
R = 20.9 Centaur 1995 SN55 and the scattered R = 23 
object 1999 RUsu- 

Of the 33 objects with well-known a in our Excited 
class, all have a > 33, though a few are Neptune-crossing 
and might be labelled Centaurs by some authors' criteria. 

Therefore if these 72 objects with good orbits are a 
guide, a few percent of all objects would be classified 
differently if full orbital elements were used instead of 
just d and i. 

3.2. Statistical Methods 

We wish to ask what forms for the differential surface 
density S(i?) = dN/dRdVl are most consistent with the 
collected survey data. Note that throughout this paper 
we will consider the differential distribution with mag- 
nitude instead of the cumulative distribution that is fit 
in most previous works. The expected number of detec- 
tions from a perfect survey over solid angle in a small 
magnitude interval Ai? is 

AN = j:{R)nAR. (4) 

Appendix IXI is a detailed explanation of the form of the 
likelihood L of observing TNOs at a set of magnitudes 
{rrii} given an assumed E(i?). This is general complex 
if the details of light curves, photometric errors, color 
corrections, and detection probability must be consid- 
ered. The Appendix demonstrates that it is safe to take 
a simplified approach that ignores many of these details, 
which we present here. 

The true surface density must be convolved with: 
the color conversion to the observed-band magnitude m; 
the measurment error on m due to noise and variability; 
the detection efficiency; and any inhomogeneities of the 
survey, leaving us with a function g{m) that describes 
the expected distribution of measured magnitudes in this 
survey. The expected number of detections from some 
particular survey is 



dmg{m)= dRnr]{R)T.{R), 



(5) 



where r/(i?) is the detection probability for a TNO of 
mean magnitude R that lies within the geometric area 
fl of the survey. This quantity can be determined from 
Monte Carlo tests. 
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Table 2. Summary of TNO Surveys 



Abbreviation 



Reference 



ioff- (deg2)- 




N (CKBO)'' 


N (Excited)*^ 


P(< N)" 




P{<C) 


0.019 


28.7 


3 





0.16 


0.65 


0.42 


0.009 


26.8 


1 


1 


0.98 


0.91 


0.09 


0.322 


25.9 


8 


9 


0.98 


0.03 


0.83 


2.30 


25.1 


17 


15 


0.49 


0.18 


0.58 


28.3 


23.8 


39 


28 


0.27 


0.44 


0.27 


296. 


20.8 


1 


5 


0.97 


0.05 


0.30 


1430. 


20.2<^ 





2 


0.28 


0.63 


0.58 



ACS This work 
CB Chiang fc Brown ^1999) 
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''Effective search area within 3° of the invariable plane at bright magnitudes, and R magnitude at which effective area drops by 50%. 
''Number of detected TNOs in the two dynamical classes defined in the text. 

•^Cumulative probabilities of this survey under the best-fit two-power-law model, for Poisson test, Anderson-Darling test, and C tests, as 
described in text. Boldface marks indications of poor fits. 
•^The TB data are not used faintward of 20.2 mag. 



The likelihood of observing a set of N magnitudes {m, } 
under an assumed distribution g{m) is 



N 



L{{mi\\g) (xe ^ Y[g{mi). 



(6) 



i=l 



In this work we will make the approximation that the 
difference between the observed magnitude rrii and the 
true R magnitude is minimal (aside from a constant color 
term) so that we may approximate 

g{m) ^ flri{'m)Y,{m) = D,cff{m)'S{'m) (7) 

N 

^ i({ma|S) a e-^ II (TO,)E(m,) (8) 



(9) 



N= dmflcs{m)'S{m). 



When fitting alternative forms of to the survey 

data, the one that maximizes the likelihood ||SJ (times 
any prior probability on the models) is the Bayesian pre- 
ferred model. Confidence intervals on the parameters of 
the underlying E(i?) can be derived from this probability 
function as well. We will always take the prior distribu- 
tions to be uniform, with the exception that the overall 
normalization of S has a logarithmic prior. 

3.2.1. Goodness of Fit 

We will be producing models for S(m), and hence 
g{m), which best fit the data. We then ask whether the 
observations are in fact consistent with having been pro- 
duced by this model. The general approach is to define 
some statistic S and ask whether the measured S is con- 
sistent with the range of S produced by realizations of 
the model. We will test goodness-of-fit with two statis- 
tics. The first is simply the likelihood L({mi}|S) itself, 
given in Equation JHl. The probability P(< L) of a re- 
alization of the model having lower likelihood than the 
measurements will be calculated by drawing random re- 
alizations from the best-fit distribution. Values P < 0.05 
or P > 0.95 are signs of poor fit. 

We als o use the Anderson-Darling (AD) statistic, de- 
fined as l)Press et al. 2002) 



AD 



[Sjm) - Pjm)]^ 
P{m)[l - P{m)] 



Here P{m) is the cumulative probability of a detection 
having magnitude <m, soO<P<l. The cumu- 
lative distribution function of the observed objects is 
S{m). The AD statistic is related to the more familiar 
Kolmogorov-Smirnov (KS) statistic, which is the maxi- 
mum of \S{m) — P(m)|, but is more sensitive to the tails 
of the distributions. We calculate Qad, the probability 
of a random realization having higher AD value than the 
real data. Values of Qad < 0.05 indicate poor compli- 
ance with the model distribution. 

Because the likelihood and AD values of the real data 
are calculated from a g{m) that is the best fit to the data, 
it is necessary to also fit each random realization before 
calculating L or AD. Because the normalization of g is 
always a free parameter, we fix each random realization 
of g{m) to have the same number of detections as the 
real data. 

We note further that we always sum the effective areas 
and detections of all surveys before analyzing the data, 
rather than considering the likelihood of each component 
survey. We believe this makes the fit a little more robust 
to small errors in individual surveys' detection efficiency 
estimates. In ^3.71 we examine whether each constituent 
survey is consistent with the S(P) derived from the full 
dataset. 

3.3. Single Power Law Fits 

Previous fits to the magnitude distribution of TNOs 
have assumed that the cumulative, and hence differen- 
tial, distributions fit a single power law. We attempt to 
fit the ACS and previous survey results to a differential 
distribution of the form 



S(P) = S23lO"(^-23) jgg- 



(11) 



dP{m). 



(10) 



We first fit this law to the older surveys, excluding the 
ACS and TB data. We recover a best fit of a = 0.61±0.04 
for the TNO sample, which is consistent with the previ- 
ous fits, e.g., TJL. This best- fit power law is a marginally 
acceptable fit to the data, with likelihood probability 
P(< L) = 0.92 and AD probability Qad = 0.06. 

When we include the ACS and TB data in the power- 
law fit we find the best-fit slope drops to a = 0.58 ±0.02. 
The fit is strongly excluded, however, P(< L) = 0.997 
and Qad < 0.001. The probability of detecting so few 
objects in the ACS survey under this power law is P(< 
N) < 10^^*, and the TB survey is also highly deficient. 
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A single power law extending to the ACS data is ruled 
out at very high significance, as expected. By contrast, 
there is 16% probability of finding < 3 TNOs in the ACS 
survey under the best-fit double-power-law I](m) f ii3.6f) . 

3.4. Binned Representation 

Figure |31 presents a non-parametric, binned estimate 
of the differential surface density. The survey data for 
each one-magnitude interval from 18 < i? < 29 are fit 
to the form l[ll|l . with a fixed to 0.6 and S23 free. The 
expectation E23 and the 68% Bayesi an credible regions 
are calculated as described in iiA.2.21 The middle panel 
shows the resulting expectation of S at the center of 
each bin, and the lower panel plots the E23 values, i.e., 
the deviation of each bin from a pure a — 0.6 power 
law. This plot is useful for visualizing the departures 
from power-law behavior, but we always fit models to 
the full survey data rather than the binned version. It 
is immediately apparent that the TNO surface density 
departs from a single power law at both the bright and 
faint ends of the observed range, for both the CKBO and 
Excited subsamples. 

3.5. Rolling Power Law Fits 

As a next level of complication, we consider a surface 
density with a rolling power law index: 

= S23l0["(«-23)+"'(i^-23)^]^ (12) 

Note that this is a log-normal distribution in the flux, 
roughly so for diameter as well. 

This fit to the full TNO sample is now acceptable, 
with Qad = 0.55 and P(< L) = 0.18. The CKBO 
and Excited samples are only marginally well fit, with 
Qad = 0.04 and 0.06, respectively. Best-fit parameters 
are given in Table |31 and the best fit are plotted 

over the data in Figure^ 

The addition of the single parameter a' to the single- 
power-law fit leads to highly significant improvements in 
the likelihood: logL is increased by 32, 22.2, and 12.6 
for the TNO, CKBO, and Excited samples, respectively. 
This is equivalent to Ax^ = 2A(logL) > 25 for one 
additional parameter, which has negligible probability of 
occurring by chance. Hence the single-power-law fits are 
strongly excluded. 

Using the Bayesian approach of H3.2l we may produce a 
probability function P{Yi23, ci, a'\{mi}) given the obser- 
vations. Figure [51 plots the credible regions for a and a' 
in the three samples. Note first that a' = is strongly 
excluded, i.e. a rolling index is required. For the rolling- 
index model, any a' < gives convergent integrals for 
TNO number and mass at both bright and faint ends. 
Second we see that the CKBO sample requires a larger 
a, meaning that is steeper at i? = 23 than for 

the Excited sample, i.e. the CKBO sample is shifted 
to fainter magnitudes relative to the Excited sample, by 
about 1 mag. In the next section we will discuss the 
implications of this magnitude shift. 

3.6. Double Power Law Fits 

We next consider a surface density that is the harmonic 
mean of two power laws: 




o 

CD 

CD 
O 
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S(i?) = (l+C)E23 



J0("2-"i)(^oq-23)^ 



(13) 

(14) 



Fig. 3. — The top panel shows the total effective survey 
area (left axis) within ±3° invariable latitude vs magnitude, 
both summed (solid curve) and for individual surveys (dashed 
curves). The histogram shows the number of detected TNOs 
for the combined surveys (right axis). The middle panel plots 
the binned estimate (Bayesian expectation and 68% credible 
range) of the differential TNO surface density near the in- 
variable plane. Solid triangles are for the full TNO sample, 
open squares (red) are for the CKBO sample, and open stars 
(green) are the Excited sample. The latter two are slightly 
displaced horizontally for clarity. The dashed line is the best 
single-power-law fit to the older data. The lower panel shows 
the binned surface density relative to an a = 0.6 power law 
(same symbols), i.e., the E23 values from a stepwise fit to 
Equation IIH . The departure of all samples from a simple 
power law is clear. This plot is also useful in that the ver- 
tical axis is the mass per magnitude interval, if the albedo, 
material density, and distance are independent of magnitude. 



Under the convention a2 < ai, the asymptotic behav- 
ior of this function is a power law of indices ai at the 
bright end and a2 at the faint end, with the two power 
laws contributing equally at i?eq- The free parameters 
for this model are {ai, a2, J?eq, ^23}- We introduce the 
double-power-law model for two reasons: first, in the 
next section we will be interested in how strongly the 
parameterization of S(-R) affects our conclusions, so we 
want some alternative to the rolling-index model. Sec- 
ond, some models for accretion/erosion of planetesimals 
predict asymptotic power-law behavior, which is absent 
in the rolling-index model. 

The double-powcr-law model adequately describes the 
TNO, CKBO, and Excited samples, with Qad > 0.12 
and P{< L) > 0.16. The best-fit parameters are listed 
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Table 3. Best Fit Differential Surface Density Models 



Sample 






Rolling Power Law 








Double Power Law 








a 


a' 






Qad 


ai 








^{< L) 


Qad 








(mag~^ deg~^) 












{mag~^ deg~^) 






TNO 


0.66 


-0.05 


1.07 


0.18 


0.55 


0.88 


0.32 


23.6 


1.08 


0.16 


0.12 


CKBO 


0.75 


-0.07 


0.53 


0.54 


0.04 


1.36 


0.38 


22.8 


0.68 


0.71 


0.23 


Excited 


0.60 


-0.05 


0.52 


0.04 


0.06 


0.66 


-0.50 


26.0 


0.39 


0.24 


0.13 



1 - ir 




R Magnitude 

Fig. 4. — Best-fit models for tlie differential surface den- 
sity of TNOs are plotted along with the binned repre- 
sentation of the data from Figure |21 The (green) stars are 
binned data for the Excited sample and (red) squares are the 
CKBO sample, and we have again divided out the function 
■|^g0.6(i?-23) ^ The two heavy curves are double-power-law fits, 
and the two thin curves are rolling-index power laws. Red 
(solid) curves are for the CKBO sample and green (dashed) 
for the Excited sample. The dash-dot curve is a secondary 
double-power law fit to the Excited which is consistent with 
these data, but inconsistent with the existence of Quaoar- 
sized objects (or Pluto). The precipitous drop in the best-fit 
double power laws at the bright (CKBO) and faint (Excited) 
ends is an artifact of the absence of detections in this sample. 
There are less precipitous drop-offs that are quite consistent 
with the data, as is apparent from Figure El 



in Table 121 and plotted with the binned representation of 
the data in Figure^ 

The values of logL for the double-power- law fits are 
within ±1.3 of those for the rolling-power-law fits. So 
while the double power law is clearly superior to the sin- 
gle power law, the likelihood itself offers no preference 
over the rolling power law. The Anderson-Darling statis- 
tic is, however, more acceptable for the double than for 
the rolling power law fits to the CKBO and Excited sam- 
ples (Table OJ- There is weak statistical preference and 
theoretical prejudice for the double power laws; in H4.2I 
we note that the rolling-power-law fits do not properly 
describe the number of very bright Excited TNOs found 
away from the invariable plane. 

In Figure we plot the Bayesian posterior distribu- 
tion P(ai,a2) for the double-power-law fits to the vari- 
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Slope a at R=23 

Fig. 5. — Allowed ranges of the slope and derivative for 
rolling-power-law fits to the sky density of TNOs, as per 
Equation dJ. Shaded regions are for the full TNO sample, 
while solid contours are for the CKBO and Excited subsam- 
ples. In all cases contours enclose 68% and 95% of the total 
posterior density. The curvature ct' of S(i?) is clearly non- 
zero, and the two dynamical subsamples are distinct. The 
lower a for the Excited sample implies that its mean magni- 
tude and mass are larger than those of the CKBOs. 



ous samples after marginalization over the less interesting 
variables S23 and i?eq- We also plot the projections onto 
the single variables ai and 02 for the CKBO and Excited 
subsamples. Note that we have applied a prior restric- 
tion —0.5 < ai < 1.5 as we consider the more extreme 
slopes to be unphysical. 

Several features of the ai, a2 constraints are notewor- 
thy. First, the CKBO and Excited samples once again 
appear to be distinct, except that the 95% CL region for 
the Excited sample has a tail at ai > 0.8, a2 ~ 0.4, that 
contains 10% of the posterior density, and overlaps the 
CKBO region of viability. Fits in this secondary range 
predict very few Excited TNOs at R < 19, which is con- 
sistent with our limited sample, but inconsistent with the 
membership of Pluto, Quaoar, and/or 2004 DW in the 
Excited class. If we include in our likelihood function 
a prior equal to the probability that each model pro- 
duce at least 1 Excited TNO at i? < 18.5 (the "Quaoar 
prior" ) , then this long tail disappears from the Excited 
credible region (as illustrated by the dash-dot lines in 
Figure El . We further discuss the CKBO/Excited di- 
chotomy in t l3.8l 

For any a2 < 0.6, the mass integral converges at the 
faint end ( ^4.1|l . and this is satisfied at high confidence 
for all samples. The faint-end slope of the CKBOs is well 
constrained at 0.38±0.12 (95% CL). The faint-end slope 
of the Excited class is poorly determined, with only a 
bound a2 ^ 0.36 (95% CL with the Quaoar prior). The 
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Fig. 6. — Allowed ranges of the two slopes for double-power- 
law fits to the differential surface density of TNOs, as per 
Equation 1131 . Shaded regions enclose 68% of posterior prob- 
ability for the CKBO and Excited subsamples (red and green 
tint respectively), with outer solid contours bounding 95% 
regions. The dashed contours are for the full TNO sample. 
Along the horizontal and vertical axes are the projected 1- 
dimensional distributions of each slope. The two dynamical 
classes have distinct magnitude distributions, with the excep- 
tion of the high-ai tail on the outer Excited contour. If we 
include a prior constraint that the Excited class contain one 
object on the sky with R < 18.5 ("Quaoar prior"), we obtain 
the dash-dot contours instead. The bright-end slope of the 
CKBO group is likely steeper than the Excited class, and the 
faint-end slope of the Excited class is probably shallower. 



absence of Excited TNOs in the ACS survey leads to this 
degeneracy. 

For ai > 0.6 the bright-end mass converges, and this 
is satisfied at 95% confidence for the Excited subclass 
and at very high confidence for the CKBO sample. The 
bright-end slope for the Excited class is 0.66]';° (95% 
CL vifith Quaoar prior) , while the absence of bright CK- 
BOs leads only to a bound of ai > 0.97 for their asymp- 
totic index. 

The CKBOs certainly seem to have a steeper bright- 
end slope (fewer large objects) than the Excited objects, 
and there is less secure evidence that the Excited class 
has a shallower faint-end slope (fewer small objects). 
There is thus evidence for different accretion and erosion 
histories for these two samples. 

3.7. Internal Consistency 

Before proceeding with further interpretation, we 
pause to ask whether there are any internal inconsisten- 
cies among the collected survey data. We check the sur- 
veys individually for consistency with the best-fit double 
power law. 

We will make use of three consistency tests. The first 
is simply the number N of detected objects. This proba- 



bility distribution for N follows the Poisson distribution. 



P{N\N) 



-N 



(15) 



and the cumulative probability P(< N) of having de- 
tected N or fewer objects is also easily calculated. The 
drawback of the Poisson test is that it makes no use of 
the distribution of magnitudes within a survey. 

The second statistic that we use is Qad statistic de- 
scribed above, which has the disadvantage that it dis- 
cards information on the total number of detected TNOs. 

The third statistic we will use is a form of log- 
likelihood: 



C{{m^}\g) = N log N -log N\- N 



N 



logff; 



J g{m) dm 



(16) 
(17) 



This statistic is useful in that it is the log of L in Equa- 
tion © when the model and number of detections N 
are fixed. The expectation value of £ when N is held 
fixed is also equal to the log of the Poisson probability 
in Equation (|15l) . This statistic hence has sensitivity to 
both the number and distribution of detections in a sur- 
vey under test. For any given model and survey, we can 
generate 1000 or more Monte-Carlo realizations to calcu- 
late the probability P{< C) of the measured likelihood 
being generated by chance under the model. 

When comparing a constituent survey to the full-data 
fits, we do not re-fit each realization, because individual 
surveys do not heavily infiuence the overall fit. 

The results of the three statistical tests for each survey 
are given in Table|3 The only sources of tension, marked 
in boldface in the Table, are for the Larsen et al. survey, 
which contains too many objects at 97% CL, and the 
Gladman et al. survey, which is overabundant at 98% 
CL. The AD tests also indicates that the Larsen et al. 
survey is too skewed toward faint objects {Qad — 0.05) 
and the Gladman et al. detections are also too skewed 
toward faint objects {Qad — 0.03). Chiang & Brown 
were also slightly lucky to find 2 objects in their 0.01 
deg^ survey, if the collective fit is correct. 

These excursions are worse than we would expect from 
Poisson statistics, but not horribly so: with 3 statistical 
tests for each of 7 surveys, we expect « 1 to show a dis- 
crepancy at > 95% significance, while we have two sur- 
veys discrepant at this level. There is no justification for 
excluding any particular survey data. For example, con- 
sider the fact that the ABM and Gladman et al. surveys 
are in poor agreement in the 25 < i? < 26 magnitude 
range. The Gl sky density in this bin is 3 times that of 
ABM. The odds of obtaining by chance a disparity this 
large given the number of detections in the survey are ap- 
proximately 1%. Since we have overlap between different 
surveys in several of our bins, the chance of our having 
found one such discrepancy between any two surveys in 
any of our bins is perhaps 5%. The discrepancy is hence 
worrisome but not outrageous. 

It is possible that either ABM overestimate their com- 
pleteness or Gladman et al. have some false-positive de- 
tections at their faint end. There are no obvious flaws 
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to either work — ABM have a thorough artificial-object 
estimate of efficiency, and Gladman et al. detect each 
object on two consecutive nights. In the absence of any 
reason to reject either dataset, we will continue to sum 
the effective areas and total detections of both surveys. 
We have verified that none of our conclusions are signif- 
icantly affected by omission of either dataset. 

Further data in this magnitude range are clearly desir- 
able, as it helps define the departure of the faint end from 
a power-law slope. Surveys at 25 < i? < 26 require long 
integrations on large telescopes with large-area CCD mo- 
saics. Such efforts are underway, using for example the 
VLT (O. Hainaut, private communication) and Subaru 
(D. Kinoshita, private communication) telescopes. 

3.8. Dynamical Subclasses 

The parametric S(i?) fits to the CKBO and Excited 
subclasses appear to differ, though the evidence is not 
yet ironclad. Since the CKBO and Excited subclasses 
have the same effective area at a given magnitude, we 
may apply a 2-sample Anderson-Darling test to see if 
their magnitudes are drawn from the same distribution. 
We obtain Qad — 0.039, rejecting at 96% confidence the 
hypothesis that the CKBO and Excited TNOs have iden- 
tical magnitude distributions. The largest difference in 
the magnitude distributions is the lack of bright CKBO 
members, which is noted by TB and discussed in detail 
below. Nearly all the statistical significance of the re- 
sult arises from the TJL sample. The test indicates that 
there are (at least) two distinct size distributions in the 
TNO population, and hence the magnitude distribution 
should be fit by dynamical class rather than summed. 

Our division into dynamical classes is crude because of 
incomplete orbital elements f H3.1|l . which can only have 
ameliorated the distinction between the two size distri- 
butions. A more precise division may yield even more 
pronounced size differences between dynamical classes. 

4. INTERPRETATION 

4.1. The Mass Budget of the Kuiper Belt 

The detection of departures from a single power law 
now make it possible to estimate the total TNO mass 
without any divergences. The total mass of a TNO pop- 
ulation may be expressed as 



TNOs 



M,, 



(18) 



= M23f7 / dR^{R)m 



-3/2 



-0.6(i?.-23) 



.0.04 



42AUy V 1000 kg m 



23 = 



:7.8 X 10^* kg 




The surface density S is the mean over solid angle fl 
of the sky, and / is the fraction of the TNO sample at 
magnitude R that lies within the area fl. The material 
density, albedo, and heliocentric distance are p, p, and d, 
and M23 is the mass of a TNO that has i? = 23 with the 
given canonical albedo, density, and distance. We ignore 
the effects of illumination phase, heliocentric vs geocen- 
tric distance, and asphericity. The angle brackets indi- 
cate an average over the TNOs at the given magnitude. 



We make the usual bold assumption that the bracketed 
quantity is independent of apparent magnitude and can 
hence be brought outside the integral in Equation ((T^ . 
We will carry the integral from 14 < i? < 31. 

4.1.1. The Mass of the Classical Kuiper Belt 

The approximation of a common heliocentric distance 
is workable for the CKBO sample, which by definition 
ranges from 38 < d < 55 AU. Of the nearly 1000 TNOs 
detected to date, none are known to have low-inclination, 
low-eccentricity orbits with semi-major axis a > 50 AU 
or a < 38 AU. A sharp decrease in surface density be- 
yond 55 AU is apparent even after correction for selection 
effect s (lAllen^^firnatfii n fc Malhotral200lUTruiillo et a,l\ 
l2nnH iTruiillo fc Brownll200lD . It is therefore physically 
meaningful to consider our CKBO sample to represent 
a dynamical class that is largely confined to heliocentric 
distances of 42 ± 10% AU. 

For the CKBO sample, the value of the integral is 
(2.85±15%) deg"^ (95% CL) when marginalized over the 
double-power-law fits. The value when using the rolling- 
index form for S(i?) is indistinguishable, so we believe 
this to be robust to parameterization. The solid angle 
over which E has been averaged is 360° x 6° . Because the 
CKBO sample is by definition restricted to i < 5°, the 
residence fraction / is high — unity for i < 3°, and 0.83 
if the inclinations are uniformly distributed between 0° 
and 5°. We will take / = 0.9 and ignore any uncertainty 
as it will be small in comparison to that of the albedo 
and density terms. We then obtain (at 95% CL) 

AfcKB = (5.3 ± 0.9) X 10^2 kg 



V0.04/ V 



42AUy VlOOOkgm-3 j*'^"^'* 

The prefactor is now determined to much greater accu- 
racy than the scaling constants. The CKBO mass is nom- 
inally equal to just O.OIOM^ or a mere 4 times the mass 
of Pluto (8 times, if the CKBOs share Pluto's density of 
2000 kg m-3). 

TJL report an estimated mass for CKBOs with diame- 
ters in the range 100 km <D < 2000 km « 1.8 x 10^^ kg, 
under the same assumptions about density and albedo as 
made here. We would naively expect our estimate to be 
larger, not three times smaller, than the TJL estimate, 
because we now include bodies smaller than 100 km — 
though the ACS data show that these smaller bodies 
hold a minority of the mass. The discrepancy is in part 
attributable to our more restrictive definition of the clas- 
sical belt: we require z < 5° whereas TJL demand only 
41 < a < 46 AU and e < 0.25. The larger part of the 
discrepancy is due to to our conclusion that the CKBO 
population has a > 0.6 for m < 24, greatly reducing the 
mass that TJL's a ~ 0.6 places in large objects. 

Gladma n et al\ (|2001) estimate the mass of TNOs in 
the 30-50 AU range to be 0.04-0.1 Af© (for unspecified 
material density), if the size distribution turns over to the 
Dohnanyi slope for I? ~ 40 km. Our new value remains 
below this estimate even if we allow for the fact that 
the Excited TNOs contribute a similar mass density to 
the CKBOs in this distance range (next section), again 
reflecting the fact that the present data fall below the 
assumed power laws at both large and small object sizes. 
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4.1.2. The Mass in Excited TNOs 

The TNOs in our Excited sample — a mixture of reso- 
nant and high-excitation non-resonant orbits — are drawn 
from dynamical families whose radial and vertical extent 
remains quite uncertain. In particular our selection ef- 
fects for the so-called scattered-disk objects are poorly 
known. Our estimates of the total mass will therefore be 
much less secure than for the CKBO population. 

The assumption that the bracketed quantity in ll9l is in- 
dependent of magnitude is dubious for the Excited class, 
because it is likely that the brighter bins are biased to- 
ward objects near perihelion, which can be an extreme 
bias for a scattered-disk member. We proceed nonethe- 
less. The integral over the double-power-law is well 
determined: marginalizing over the fitted parameters we 
obtain (3.5 ± l)deg~^ for the double power law, some- 
what lower for the rolling- index fit. 

Because the excited population has high inclinations, 
the residence fraction / will be smaller; if the inclina- 
tions i are distributed as a Gaussian with ai = 10°, then 
they spend on average / « 0.5 of their time within 3° 
of the invariable plane. The inclination distribution is, 
however, poorly known, especially under our definition 
of the Excited sample. 

The mean distance appears as dP and is highly un- 
certain for the "scattered-disk" objects, but is sensibly 
bounded for Plutinos. An upper bound on the Plutino 
mass comes by assuming that the Excited class contains 
most of the Plutinos that are closer than 39 AU to the 
Sun, or very crudely half of the total Plutino population. 
The brighter surveys are biased against Plutinos that are 
currently beyond their semi-major axis distance, and un- 
der our definition of the CKBO sample, Plutinos at low 
inclination beyond 38 AU will be put into the CKBO 
class. So we will calculate a total mass for the detectable 
Excited sample assuming d = 39 AU, and double it to 
bound the Plutinos. Following this procedure, we obtain 
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An estimate for the so-called "scattered-disk" popula- 
tion is even less certain given the potentially large val- 
ues of (d^) (even though no detected objects^ yet exceed 
d ~ 60 AU). Taking the mean d to be 42 AU and / = 0.5 
gives an Excited-class mass of 1.3 x 10^'^ kg. It perhaps 
suffices to estimate that the Excited sample is compara- 
ble to or several times the mass of the CKBOs, because 
the surface-density integrals are nearly equal, the resi- 
dence fraction / is lower for the Excited objects, and 
the mean d^ factor could be larger for the Excited class. 
The value here is, however, a fair estimate of the Excited 
mass within 50 A U. 

It thus appears that Pluto itself accounts for > 10% 
of the mass of the Plutino population, and perhaps of 
the entire Excited population. We will investigate this 
further in H4.2I 

TJL crudely estimate the number oi D > 100 km ob- 
jects with scattered-disk orbits to be comparable to the 

^ We do not include the newly-discovered object Sedna at d = 
90 AU in the scattered disk given its perihelion at 76 AU. 



CKBOs, and the Plutinos to be ^ 20 times less abun- 
dant. The absence of Plutinos in the ACS data suggests 
that the Plutinos also have a break in S(i?) such that lit- 
tle of their collective mass is in small objects, but we do 
not attempt a further quantitative bound. If this TJL 
estimate is correct, then Pluto dominates the Plutino 
mass and is several to tens of percent of the scattered- 
disk mass. 

4.2. The Largest Objects in the Dynamical Classes 

We examine here some implications of the large-size 
behavior of our fitted S(-R) functions. These results are 
not of course a direct consequence of the new small-end 
ACS data, except insofar as the latter make it clear that 
single-power-law fits to the full population should not be 
expected to properly characterize the bright end. The 
results of this section flow from our efforts to create a 
homogeneous subsample of all the published brighter sur- 
veys. 

The double-power-law fits to the CKBO sample favor 
a bright-end slope that is significantly steeper than for 
the Excited sample. If correct, this implies a different 
accretion history for the two classes, since it is believed 
that erosion should not have affected the largest TNOs 
('Sternll995D. It would also imply that the brightest (and 
largest) member of the Excited class will be significantly 
larger than the largest CKBO, and that the Excited sam- 
ple will have a substantial fraction of its mass in the 
largest objects (since ax ~ 0.6) while the CKB will not. 
A possible preference for large objects to reside in high- 
excitation populations has been s uggested before in stud - 
ies of less well-controlled sa mples (Levison fc Sterni2C)0ll) 
and in theoretical studies l)Gomes.2003.'l . and is verified 
by the TB survey. 

We have inferred the bright-end slope difference by us- 
ing only the subset of survey data to date that has invari- 
able latitude < 3° and has published survey statistics. 
We may check the accuracy of our inferences by some 
comparisons with the full sample of nearly 1000 TNOs 
that have been discovered to date, as listed on the Minor 
Planet Center web pages July 2003. 

We first ask what should be the brightest TNO on the 
full sky if extrapolation of our fitted double-power-law 
I](i?) were to describe correctly all TNOs. We may cal- 
culate the probability distribution for the magnitude of 
the brightest TNO by marginalizing over all our double- 
power-law fits, and assuming an available area of sky 
equal to ^/ f, using the values estimated above. To be 
specific, we calculate 

P[N{< R) = 0]=J dp P(p|{mJ) P[Ni< R) = 0\m) 
dpP{p\{m,}) 



exp 



^{n/f) / dm'l](m'|p) 



(24) 



Here p is a vector of parameters for S, P{p\{mi}) is the 
normalized Bayesian posterior probability for the model 
parameters given the observed data, and = 360° x 6° is 
the total area of the low-invariable-latitude strip on the 
sky. 

Figure [7| shows our estimated probability of there be- 
ing, somewhere on the sky, a TNO brighter than a given 
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R magnitude 

Fig. 7. — The distribution for the brightest expected TNO 
on the sky is shown vs R magnitude for the CKBO (red, 
right) and Excited (green, left) samples. The plotted quantity 
is the cumulative likelihood of finding a TNO brighter than 
R. The solid lines are derived from an extrapolation of the 
double-power-law fits to the invariable-plane surveys collected 
in this paper, and the dashed lines are for the rolling-index 
power laws. Given that much of the relevant sky area has 
now been searched for bright TNOs, we see that the extrapo- 
lations of the power laws are consistent with 2000 KX14 being 
the brightest CKBO (or nearly so) , and with Pluto or Quaoar 
being the brightest Excited TNO [Pluto is plotted at the ap- 
parent magnitudes it would have with albedo 0.04-0.10.] If 
the rolling-index fit is the proper extrapolation, Pluto would 
have to be considered an exceptional object. 



R magnitude. For the CKBOs, the double- and rolling- 
power-law fits are in general agreement. The median ex- 
pected magnitude for the brightest CKBO on the sky is 
20.3 mag. The brightest TNO found to date that meets 
our CKBO criteria is 2002 KX14, at i? = 20.6, which is 
found by the TB survey but was not used in our analysis 
since it is below our adopted 20.2 mag cutoff. Nearly all 
the sky within 5° of the invariable plane has been sur- 
veyed for such bright objects, by TB and/or the Southern 
Edgeworth-Kuiper Survey (Moody et al. 2003), with all 
discoveries having been transmitted to the MFC. Hence 
it is unlikely that a CKBO significantly brighter or larger 
than 2002 KX14 exists on the sky. The fitted and ex- 
trapolated S(-R) models suggest a 70% chance of finding 
a CKBO brighter than 2002 KX14, so it would be ac- 
ceptable for 2002 KX14 truly to be the largest CKBO 
(or nearly so). 

For the Excited class, an extrapolation is more spec- 
ulative since a smaller fraction of the available sky has 
been surveyed and we include in our analysis only a frac- 
tion of the surveyed area. The median expected magni- 
tude of the brightest Excited TNO is 18.4 for the dou- 
ble power law, and 19.0 for the rolling-index fits. The 
brightest TNO is Pluto, which would h ave R = 17.5 if 
it had our assumed albedo of 0.04 (Bur atti et oi.ll2(ffl^ . 
OT R = 16. 5 if the albedo is 0.10, comparable to Quaoar 
ijBrown fc Truiillo 2004) . Under the double-power-law 
fits, there is 4-17% chance of finding an Excited TNO 



this bright, so it is feasible to consider Fluto to have 
been produced by the same physical process as the other 
TNOs. Under the rolling-index model, however, the 
odds of a Fluto on the sky are < 2%, hence under this 
model we would have to consider Fluto as an excep- 
tional object with a distinct formation mechanism from 
the other TNOs. The largest and brightest TNO out- 
side the Fluto/Charon system is Quaoar at i? = 18.5 — 
also discovered by TB but outside the bounds of our 
±3° swath — which at ? = 8° would be excluded from our 
CKBO sample. From the Figure we see that if the bright 
end follows a power law, we would expect the largest 
Excited TNO to be in the Quaoar-Fluto range. So we 
again find that the largest objects known to date are con- 
sistent with extrapolation of our double-power-law S(i?) 
fits — or with the rolling-index model, if we consider the 
Fluto-Charon system as exceptional. 

The degree of dominance of the brightest known ob- 
jects is also consistent with our double-power-law fits to 
our selected subsample. The ten brightest known CK- 
BOs have 20.6 < R < 21.8 and contribute « 2% of the 
estimated CKB mass (assuming common albedo and ma- 
terial density). Again, surveys for low-inclination TNOs 
this bright are now majority complete, so this fraction is 
not likely to evolve much. Hence the largest objects in 
the CKB do not hold a significant fraction of its mass, 
confirming the validity of the ai > 0.85 bright-end slope. 

The ten brightest known resonant/scattered TNOs, ex- 
cluding Fluto and Charon, have 18.5 < R S 20.5. Their 
combined mass, with canonical albedo and density, is 
1.3 X 10^^ kg, roughly equal to Fluto's mass, and roughly 
10% of the total estimated Excited family mass. Thus 
the largest bodies known hold tens of percent of the 
known Excited-class mass, confirming that the power- 
law slope is near 0.6. 

The conclusions we have drawn from our controlled 
TNO subsample are therefore supported by analysis of 
total sample of known bright TNOs. In particular we 
find that the Excited class is near the ai = 0.6 value 
with equal mass per logarithmic size bin, while the CK- 
BOs have steeper ai > 0.85 that puts less mass in 
large objects. Under the doublc-power-law model, Fluto 
is a uniquely but not anomalously large Excited TNO. 
Since erosion should not have significantly altered the 
^ 100 km objects in the current trans-Neptunian envi- 
ronment, a simplistic interpretation is that the Excited 
objects were formed in a region of the Solar System that 
allowed the accretion process to proceed further than in 
the ^ 42 AU location of the CKBs. This is consistent 
with the scenario that the Excited TNOs formed in a 
low-excitat i on population closer to the Sun (Malhotr^ 
'1993, 1995; Gomes 2003), where higher space densities 
and relative velocities could speed accretion. 

4.3. The Smallest Known TNOs and Erosion History 

Calculations of collision rates and energies suggest that 
the population of TNOs with D < 100 km has been 
strongly influenced by coUisional erosion. The ACS sur- 
vey provides the first strong constraints on this popula- 
tion, so in this section we provide some basic comparison 
of the observed magnitude distribution to models of the 
accretion/erosion history of the TNOs. Throughout this 
section we will make the simplifying assumption that the 
TNOs samples we have gathered are at fixed distance 
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with fixed albedo, so that the magnitude distribution 
may be directly translated into a size distribution. If the 
magnitude distribution scales locally as dN/dra oc 10"™, 
then the size distributions will locally follow the com- 
mon parameterization dN/dD cx D^'' with q = 5a + 1. 
We will prefer the double-power-law models for to 
the rolling-index models for this discussion, as the size- 
evolution models typically predict such behavior. 

4.3.1. Comparison with the MFC Database 

The conclusion that the faint-end is shallower 

for CKBs than Excited objects does not have exception- 
ally strong statistical support in our sample, e.g. 3 vs 
detections in the ACS survey. We again look to see 
if the trend is borne out by the sample of all objects 
submitted to the MPC. Of the 23 objects with H > 9.6 
and perihelion q > 28 AU, only 5 would fall into our 
Excited class. The dominance of CKBOs in this range 
is significant given that: Excited TNOs dominate the 
\ow-H objects are comparable to CKBOs near R « 23; 
discovery of high-i/ Excited objects is favored because of 
their nearer perihelia; and several of the H > 9.6 Excited 
TNOs could be CKBOs once more accurate inclinations 
are obtained. This test is not independent of our more 
controlled analysis because many of the H > 9.6 TNOs 
at the MPC are from the ACS, ABM, and Gl surveys. 
Nonetheless the trend in our controlled subsample seems 
to extend to the full sample. 

4.3.2. Equilibrium Collisional Cascades 

An analytic description of a collisional cascade is given 
bv IDohnanvil ljl969f) . who derives q = 3.52, a = 0.50, 
for an equilibrium solution. In such a steady state, the 
loss of objects from a given mass bin due to grinding 
and catastrophic destruction is exactly replenished by 
the gain of fragments from collisions on larger bodies. A 
population of high-mass objects is of course required as 
a mass reservoir to maintain such an equilibrium. We 
see from Figure that a2 > 0.5 has a low probability 
of describing either the CKBO or Excited populations 
in the ACS survey: the true distribution is markedly 
shallower (meaning relatively fewer small bodies) than 
the canonical Dohnanyi value, particularly in the Excited 
sample. 

The failure of the analytical model could perhaps be 
ascribed to the particular assumptions in the model, e.g., 
that the size distribution will follow a power law and that 
the distribution of fragment sizes will take a particular 
power-law form, although the equilibrium form should 
be robust a s long as the fragmentation law is scale-free. 
IPan fc Saril l)2004j) suggest that the observed size dis- 
tribution can be explained by a fragmentation law for 
rubble-pile objects which is not scale free. An alterna- 
tive explanation is that the present-day TNO population 
is not in collisional steady state. Perhaps we have in- 
stead a snapshot of two populations for which the small- 
body population is continuously decreasing and erosive 
destruction is very advanced. In the CKB, it is not even 
clear what bodies could serve as a mass reservoir, since 
the largest bodies hold little of the mass. 

4.3.3. Numerical Models: Ceneral Results 



Modelling of non-equlibrium collisional evolution re- 
quires numerical simulation, particularly if coeval ac- 
cretion is modelle d as well. T hree groups have pro- 
duced such models: iDurda fc Ster n (2000) (who estimate 
body lifetimes but do not follow the long- term evolu- 
tion of the size distribution') . IKenvon fc LuJ lfl999') , and 
iDavis. Farinella. Sz. Weidenschillin3 ~ ilT99^K and th e an- 
tecedents of these papers, particular Iv ISternI l)1995l) and 
IDavis fc Farinella (1997), offer results of numerical sim- 
ulations o f TNO siz e distributions. A review is given by 
iFarinella. Davis, fc St ern (2000). 

Before comparing the numerical models to the new 
data in detail, we first review some of the general con- 
clusions of these works, and re-examine them in light of 
the new observational results, namely: the detection of a 
strong break to shallower distributions a,t D < 100 km, 
implying much lower impact rates; the stronger bounds 
on the bright-end a; and the clear detection of size- 
distribution differences between two dynamical popula- 
tions. 

A simple analytical result of iSternI l)1995(l is that for 
a population with mean eccentricity 0.03 < (e) < 0.1 
at ~ 42 AU, collisions are on average erosive for objects 
smaller than a critical diameter D* and accretional for 
larger bodies. The critical diameter D* is 100-300 km 
for the "strong" bodies of Durda & Stern (2000), which 
appear at R magnitudes 23-25.5 under the canonical 
albedo. Note from Figure 0] that this is the magnitude 
range where T,{R) develops a significant deficit relative to 
the bright-end power law. The TNOs that are "missing" 
are those that are suscept ible to collisi onal destruction. 

A general conclusion of lSternI jl995fl and later models 
is that the collision rate in the present trans-Neptunian 
region is too low for D > 200 km objects to have been 
formed by pairwise accretion in 4 Gyr. The ACS results 
only exacerbate this difficulty since the mass in small 
bodies is found to be much lower than all the models had 
assumed. We note, however, that the maximum size of 
bodies in the CKB sample is significantly lower than that 
in the Excited sample, with the largest known CKBO 
being 60 times less massive than Pluto. These largest 
bodies should be largely unaffected by erosion according 
to all models. We can conclude that the Excited popula- 
tion is "older" than the CKB in the sense of cumulative 
number of accreting collisions. This is consistent with an 
origin for the Excited bodies a smaller heliocentric dis- 
tances, in a denser and faster-moving section of the disk 
t han the CKB. 

IDurda fc Ster n (2000) calculate impact and disruption 
rates for TNOs in the present Kuiper Belt, concluding 
that all TNOs should be heavily cratered but that the 
lifetime against catastrophic destruction for D > 1 km 
bodies exceeds the age of the Solar System. We now be- 
lieve the space density of I? < 100 km bodies to be much 
lower than assumed in these models, perhaps by several 
orders of magnitude for km-scale bodies. A recalculation 
of present-day impact rates could show that small (km- 
scale) TNOs have negligible cratering in the past billion 
years, and even large TNOs would have recent craters 
only from sub-km impactors. This would argue against 
a collisional resurfacing effect as the source of color di- 
versity among TNOs. 

So we find the apparent paradox that objects small 
enough to be subject to collisional disruption are strongly 



19 



depleted, yet have lifetimes against disruption that are 
longer than the age of the Solar System. A possible reso- 
lution of this paradox is that we now see the end state of 
the erosion process: the original CKB (and perhaps the 
precursor region of the Excited population) was more 
massive and richer in small bodies, leading to strong 
erosion and depletion. The erosion lowered the volume 
density to the point where lifetimes exceeded 10^ years, 
leaving the present "frozen" population behind. In this 
scenario, present-day D < 100 km bodies are likely all 
collision fragments, but they may have had little alter- 
ation in the past 3 billion years. Surfaces would be heav- 
ily cratered but old. 

At some point in the history of the CKB, the Excited 
population started to cross the CKB, with higher colli- 
sion velocities, which would have exacerbated the erosive 
depletion process. 

Alternatively, the accretion process in the original 
source disk was extremely efficient at collecting 1-10 km 
bodies into D > 100 km bodies — but not at producing 
1000 km bodies from these. A quantitative explanation 
for the preferred scale of accretion would be required. 

4.3.4. Comparison with Detailed Models 

iKenvon & LuJ ||1999D and 

IPavis. FarincUgu fc Weidenschilhnel l) 19991) plot the 
time history of the TNO size distribution in their 
numerical accretion/erosion models. The following seem 
to be generic features of the models: at the large end, 
the distribution becomes shallower with time as the 
size of the largest accreted objects increases. At the 
small end, the slope either approaches the Dohnanyi 
value iKenvon fc Luul I1999D or beco mes progressively 
shallower ijDavis.^^mnena. &: Weiden schillina 1999). 

The observed CKB and Excited distributions disagree 
with the Kcnvon & Luu (1999) models in several re- 
spects: first, the observed faint-end slope does not match 
the Dohnanyi value they predict; the models predict a 
large excess of bodies at the transition region, which is 
not observed; and the observed transition appears to be 
at the 10-100 km scale, not the km scale predicted by 
the models. These models, however, did not attempt 
to model the effect of several Gyr of further coUisional 
erosion after the formation epoch. 

The a greement of the da ta with the 
iDavis. Far inella. & Wcidcnschillin^ l|1999D models 
is fairly strong. The last time slice (1 Gyr) of their simu- 
lation shows bright- and faint-end slopes of a « 0.6 and 
« 0.0, respectively, with a transition near « 50 km. 
This is remarkably similar to the observed Excited 
sample magnitude distribution. 

The CKB sample resembles younger time slices in these 
models (~ 200 Myr) in having steeper slopes at the 
bright end (ai « 1.2) and faint end (a2 ~ 0.4), as well 
as a smaller maximum object size. The transition re- 
gion is at smaller diameters in the models than in the 
data, however. It appears in any case that the CKB is 
less evolved than the Excited sample in an erosion sense 
as well as an accretion sense. This might be expected 
from the smaller velocity dispersion of the present-day 
populations (if erosion postdates any migration of the 
Excited bodies) or from the higher space density and/or 
velocity of the Excited source region (if erosion predates 
migration). 



The ACS data should serve as a target for future size 
evolution models. It may be particularly important for 
future models to consider the coupled evolution of the 
two (or more) co-spatial dynamical populations that ap- 
pear to exist currently and to have different size distribu- 
tions. An interesting question is whether any dynamical 
mechanism is required to reduce the mass of the CKB 
to its present ~ O.OIM® value from the ~ lOM® levels 
that appear necessary to support large-object accretion, 
or whether erosive removal is sufficient. This question is 
particularly interesting in light of the absence of detected 
TNOs beyond 55 AU, as discussed below. 

4.4. The Source of the Jupiter- Family Comets 

It is currently widely accepted that the Jupiter fam- 
ily short period comets (JFCs) and the Centaurs are 
objects that have escaped from the Kuiper Belt. Our 
HST-ACS survey provides an estimate of the popula- 
tion of cometary precursors in the Kuiper Belt, which we 
can now compare with estimates obtained from dynam- 
ical models of the Kuiper Belt-JFC connection. In one 
of the earli est modern models of the Kuiper Belt-JFC 
connection, iHolman fc WisdomI l|199,'^ estimate that a 
population of 5 x 10^ cometary precursors in the Kuiper 
Belt in the 30-50 AU range is required to account for 
the observed population of JFCs. A more detailed 
calculation of essentially the same physical model by 
iLevison fc DuncanI l)1997j) revises this estimate to 7x 10^. 
These calculations are based on what can now be de- 
scribed as a "cold" Kuiper belt, with most KBOs in low 
eccentricity, low inclination orbits. In another model, 
iDuncan fc Levis on (1997) postulate the "scattered disk" 
as a source of JFCs, and estimate a population of 6 x 10® 
cometary precursors in the entire scattered disk, of which 
~ 1.4 X 10® are at heliocentric distances 30-50 AU. In yet 
another calculation, MorbidcUi ( 1997) models the Pluti- 
nos as a source of the JFCs, and estimates a required 
source population of 4.5 x 10® Plutino comets. A com- 
monality of these models are the assumptions that the 
various classes of the Kuiper Belt constitute stable reser- 
voirs of cometary precursors, and that some fraction of 
these objects escape into the inner solar system on gi- 
gayear timescales due to slow orbital chaos induced by 
the gravitational perturbations of the giant planets. The 
models differ only in the initial conditions of the Kuiper 
Belt comets, i.e., in the choice of the KB subclass for the 
putative source of the JFCs. 

We can convert the modeled population estimates to a 
surface density by assuming that the projected sky area 
of these estimated populations is 10^ deg^ (corresponding 
to a ±15° latitudinal band around the ecliptic or invari- 
able plane). A greater challenge is to define the mag- 
nitude range of cometary precursors that correspond to 
the dynamical models. CLSD took this to be i? ~ 28.5, 
corresponding to an object of ~ 10 km radius at helio- 
centric distance 40 AU, assuming the usual 4% geometric 
albedo. With this assumption, CLSD found agreement 
between their measured surface density and the required 
cold precursor population. 

We find a surface density three orders of magnitude 
lower than the claim of CLSD. The nuclei of JFCs have 
now, however, been measure d to have d iameters predom- 
inantly in the range 1-10 km l)Lamv et alJZOOA) . It is un- 
fortunately unclear whether the JFC population is close 
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Fig. 8. — The cumulative surface density of the Excited 
KBO population (green) and of the classical population (red), 
parametrized by the double power law (Equation The 
regions are bounded by the 95% confidence upper and lower 
bounds at each magnitude. Dot-dash line is the now-outdated 
power-law fit to the total population from TJL. The horizon- 
tal bands in the upper right are theoretical estimates based 
on models of the Kuiper belt as a source of the JFCs; the 
assumed source population for each case is labelled. The hor- 
izontal extent is a reminder that the precursor population 
of known JFCs is not currently well defined, probably lying 
somewhere above 1 km diameter. The observed populations 
fall short of the surface densities calculated for putative JFC 
reservoirs unless a favorable extrapolation of the Excited-class 
density is assumed. The CKB population falls ~ lOx short 
of the required reservoir even with favorable extrapolation. 



to complete for D > 1 km. It is unlikely that the JFC 
precursors would be smaller than the JFCs themselves, 
or that albedos could be significantly lower than 0.04, 
so we may consider i? « 35 to be the faintest possible 
precursor population. The dynamical estimates of the 
surface density of trans-Neptunian cometary precursors 
are shown in Figure |S| by the horizontal bands in the 
upper right, which indicate a range of 1-10 km as the re- 
quired size of the true precursors. Figure|Slalso plots the 
95% likelihood of the parameterized fits to observations, 
as described in 

Extrapolation of the ACS measurements to fainter 
magnitudes is required to compare the TNG densities 
to the putative required JFC source populations. The 
classical belt and Plutino populations are insufficient to 
supply the JFCs even under the most favorable extrap- 
olations of the ACS observations. Of course we cannot 
rule out an upturn in the T,{R) slope at i? > 29. A slope 
near q = 3.5 (a = 0.5) for 1-20 km diameter Jupiter- 
crossing bodies is suggested by Zahnle et al. ( 2003) based 
upon crater counts on Europa and the bright terrain of 
Ganymede. This value is outside the allowed a2 range 
for the Excited TNOs an d marginal fo r the C KBGs. The 
ecliptic-comet samples of lLamv et al\ ()2004|) show a size 
distribution with a = 0.32 ± 0.04 or 0.38 ± 0.06, depend- 
ing upon the inclusion of NEGs with "cometary" orbits, 
in the range I? > 3.2 km. Either range is marginally 
compatible with our Excited TNO sample, and agrees 



well with the measured a2 for CKBGs. Gf course many 
comet nuclei are known to undergo dramatic transfor- 
mation or disintegration upon entering the inner solar 
system, so the match with TNG sizes may be fortuitous. 
The apparent discrepancy between our faint TNO size 
distribution and the crater sizes on the Galilean satellites 
requires further investigation, though it is interesting to 
note that both the TNO data and the crater counts have 
size distributions that flatten to smaller sizes. 

The scattered disk remains a viable precursor popula- 
tion, if the precursors are near 1 km diameter and the 
slope at i? > 29 is the steepest permitted for i? < 29. 
Thus with a favorable choice of precursor diameter and 
favorable extrapolation to Rk, 34, present scenarios for 
JFC origin in the scattered disk can work. 

There are several factors of 2 (or more) still indeter- 
minate in the quantitative understanding of the origin 
of JFCs: the completeness of the known population, the 
size and albedos of their precursors, the dynamical de- 
livery rate and lifetime of the JFCs, the extent of the 
scattered-disk, and the extrapolation to fainter magni- 
tudes. Any definitive judgement on the origin of JFCs 
would clearly be premature at present. 

4.5. Constraints on the Distant Population 

The ACS survey fails to detect any object further than 
42 AU despite extending the TNO completeness limit 
by 2 mag from previous work. Several previous works 
llAllen. Bernstein fc: Malhotrall200l]: iTruiillo et"aIll200H 
iTruiillo fc Brownll2001l^ have shown that the space den- 
sity of objects of size D > 150 km must decrease sig- 
nificantly beyond 50 AU from the Sun, i.e., the CKB is 
bounded. It is possible that there is an outer belt that 
is deficient in objects of this size and hence undetected. 
We provide here a crude bound on an outer disk com- 
posed solely of smaller bodies based upon the complete 
absence of distant objects at i? < 28.7 in the 0.02 deg^ 
of the ACS survey. 

Take a simple model in which an outer belt is localized 
to heliocentric distance d and contains a population of 
objects of diameter D. We assume that the outer belt 
members have albedos and material densities similar to 
the CKB, and that the projected mass density on the 
invariable plane is some multiple g times the value for 
the CKB derived in Equation (EB- For d = 60 AU, 
the magnitude limit for the ACS survey corresponds to 
D > 37 km (for p = 0.04). The invariable-plane number 
density of an outer disk of such objects would be 900g per 
square degree. The ACS data limit the density at 95% 
CL to be < 150deg~^, implying g < 0.17. If D of the 
outer belt is larger, the limit on g is roughly unchanged, 
because at brighter magnitudes the area surveyed to date 
is larger. The rough conclusion to be drawn is that an 
outer belt must either have most of its mass in objects 
with D < AO km, or must be substantially less massive 
than the known CKB. Caveats are that the presumed 
outer belt must be near the invariable plane at the ACS 
field, must not be more dispersed vertically than our def- 
inition of the CKB, and must not have albedo or material 
density dramatically different from the CKB. 

This limit on the mass of an outer Kuiper Belt does 
not apply to objects with orbits like Sedna which at 
d = 90 AU is approachi ng its perihelion near 76 AU on 
a highly elliptical orbit ijBrown. Truiillo. fc Rabinowitzl 
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l2nf)l . While the ACS survey would easily detect the 
motion of bodies anywhere on such orbits, they would 
be > 500 AU distant most of the time, where even ob- 
jects 1000 km in diameter would fall below our flux We 
also do not usefully constrain the mass contained in rare 
large objects such as Sedna itself. 

A dearth of CKB mass beyond ~50 AU may be a 
result of the initial conditions of Solar System forma- 
tion. Wcidcnschillins (2003) attributes a lack of CKB 
mass to processes from the era of planetesimal formation 
4.5 billion years ago. In his model, solid bodies of sizes 
0.1-10 m drift inward from locations of 50-100 AU due 
to gas drag. Interior to 50 AU, sufficient solid surface 
densities exist for accretion to continue. Large bodies 
do not accumulate outside of ~50 AU, and the exact 
location of the outer edge depends on the initial disk 
density radial profile. This model adequately produces 
the observed CKBO distribution and implies that there 
should be little mass beyond 50 AU, even for bodies 
with D < 40 km. Alternately, the outer region of the 
nascent Solar System could have been eliminated entirely 
through photoevaporation of the circumst ellar disk due 
to UV flux from nearby OB stars {e.g., iThroop et al\ 
((2001)). Evaporation of the disk (and entrained grains) 
occurs to radii as small as 50 AU in 10^ years in an 
Orion-type environment. The 10* or 10^ year CKBO for- 
mation timescale would therefore preclude formation of 
bodies outside the e vaporation radius of perhaps 50 AU. 
iLevison fc Morbidel li ( 2004) suggest that the full TNO 
population accreted within « 30 AU and was transported 
to its present location by interaction with a migrating 
Neptune; in this case we still need to explain the trunca- 
tion of the original planetesimal system at 30 AU. 

The lack of CKBOs beyond 50 AU may instead be 
a result of dynamical evolution of the outer Solar Sys- 
tem s ubsequent to its formation. Adams & Laughlin 
l)2001t) considered the case in which the Solar System 
had dynamical encounters with other members of its 
birth cluster; these interactions could have disrupted the 
KB outside of 50 AU and increased KBO eccentricities 
to 0.2 or larger. Later close stellar encounters could 
have had similar effects: disrupting, exciting, or even re- 
moving the outer CKB l|Ida. Larwood. fc BurkerHl200?i 
ILarwood Kal^l200l. 

The observed data is also consistent with the interpre- 
tation that signiflcant mass is present in the CKB out- 
side of 50 AU and that this mass resides entirely in bod- 
ies smaller than ^-^40 km. It is conceivable that, within 
the standard model of accretion of planetesimals in a 
(gaseous) circumstellar disk, accretion timescales beyond 
50 AU were simply too long for bodies larger than 40 km 
to form before the present day or some disruptive event. 
We note that a substantial mass in D > 3 km bodies in 
an outer Kuiper Belt shoul d be detectable by o ccultation 
surveys in the near future (jAlcock et QiJl2003l) . 

A final possibility is that large bodies formed in distant 
orbits but have since been largely eroded or shattered to 
smaller sizes. It is usually difficult to destroy the largest 
bodies in a population, but the remaining ones could 
retain a sufficiently small portion of the original nebular 
mass that they remain undetected. 



The superior coverage and efficiency of the Advanced 
Camera for Surveys has allowed a survey for TNOs with 
m < 29.2 over 0.02 deg^. The faint limit of the survey is 
determined by the criterion that a TNO have signal-to- 
noise ratio ly > 8.2 over the 22 ksec total integration time 
of the discovery epoch; false positive detections begin to 
appear below this threshold. A modest investment in 
present-day CPU power allows us to track, in effect, the 
HST along each of lO" possible TNO orbits. 

Three new objects are detected, compared to the « 85 
that were expected from extrapolation of the power-law 
fit to brighter TNO survey detection rates. Artificial 
TNOs implanted into the data are recovered efficiently, 
verifying the completeness of the detection process, and 
the reality of the three objects is confirmed by their de- 
tection in a second set of ACS observations and the de- 
tection of the brightest one {R « 27) at the Keck Obser- 
vatory. 

All three detected objects are consistent with nearly- 
circular, low-inclination orbits near a ~ 42 AU, although 
complete orbital elements will require foUowup detection 
(with HST)^. Assuming an albedo of 0.04, the detected 
objects range from 25-44 km in diameter, and are all well 
above the survey detection threshold. 

No objects are detected beyond 43 AU, though we have 
the astrometric sensitivity to detect Solar System mem- 
bers many hundreds of AU distant, and photometric sen- 
sitivity to see 37 km objects at 60 AU. This further tight- 
ens the limits on any "outer" Kuiper Belt: any outer 
component must be significantly less massive than the 
known population, or predominately in the form of ob- 
jects smaller than 40 km. 

In an assemblage of data from ground-based surveys 
with well-quantified completeness in magnitude and in- 
variable latitude, a crude division of the sample into 
"classical" and "Excited" dynamical classes reveals a 
difference in magnitude distribution at 96% confidence. 
Since the two samples have roughly similar distance dis- 
tributions, we can infer a true difference in size distribu- 
tion and hence in accretion/erosion history. Both pop- 
ulations are grossly inconsistent with a singic-power-law 
magnitude or size distribution. Double power laws give 
good fits to the present observations, with rolling power 
laws slightly less acceptable. The parametric fits indicate 
that the Excited population has more of its mass in large 
objects, so that the brightest TNOs are almost entirely 
in the Excited class, as has been noted before. More 
quantitatively, we find that under the double-power laws 
which fit our limited sample, it is feasible to have a sin- 
gle body as massive as Pluto, i.e., Pluto is uniquely but 
not anomalously large under this distribution. This is a 
sign that the formation of Pluto follows the same phys- 
ical mechanism as the other Excited TNOs. The accre- 
tion history of the Excited bodies placed a larger fraction 
of the mass into the largest bodies, whereas the largest 
CKBO is « 60 times less massive than Pluto. This is 
crudely consistent with scenarios in which the Excited 
bodies were for med at smaller heliocentr ic distances than 
the CKB ( Mal hotralllQM IT995t lGmIe3 l2003) . 

The total mass of the CKB (as we have defined it, with 
i < 5°) is constrained to O.OIOM^, with accuracy now 
limited solely by estimates of the albedo, distance, and 



5. SUMMARY 



^ Cycle 13 observations on HST have been allocated for retrieval. 
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material density of the bodies. This is significantly lower 
than previous estimates, and exacerbates the difference 
between the present mass and the ~ lOM^ initial mass 
that is believed to be required to facilitate accretion of 
the larger TNOs. The total mass of the Excited sample 
is less certain because its vertical and radial extent are 
poorly known; but the total mass of objects within 50 AU 
appears similar to the CKB mass. 

The detection of features in the magnitude (and pre- 
sumably size) distribution gives fundamentally new in- 
formation to be used in constraining the history of 
the trans-Neptunian populations, particularly the ac- 
cretion/erosion history. Both dynamical samples show 
breaks to shallower size distributions for the _D < 100 km 
population. A reevaluation of the models for coUisional 
evolution is necessary; the assumption of coUisional equi- 
librium with a scale-free fracturing law may need to be 
abandoned. Estimates of cratering and resurfacing rates 
in the trans-Neptunian region will also require revision. 
It is possible, for example, that collision rates are too low 
to allow significant exposure of fresh ices, a mechanism 
that h as been invoked to e xplain the color diversity of 
TNOs l|Luu fc .TewittllTool . The CKB seems to be ex- 
cluded as a source of the JFCs, with the Excited class 
being a viable JFC reservoir only with optimistic choices 
for model parameters and TNO extrapolations. 

The present data are consistent with a picture in which 
both the CKB region near 42 AU and a now-vanished 
inner region of planetesimals were originally much more 
massive than at present. Accretion of 10 km bodies pro- 
ceeded in both regions, advancing to 100 km bodies in the 
CKB; 1000 km bodies are produced in the inner region, 
where higher surface densities promote more rapid accre- 
tion. The migration of Neptune (or some other massive 
body) clears the inner region, scattering a few percent 
of the bodies into present "Excited" orbits. This mi- 
gration or some other process also excites the eccentric- 
ities of the disk, commencing erosion among TNOs with 
D < 100 km. These smaller bodies are greatly depleted, 
today's population being only a fraction of the fragments 
produced in early collisions. Meanwhile, beyond 50 AU 
there is either little solid material, accretion never pro- 
cedes beyond 10 km bodies, or the modest-sized bodies 
are destroyed or removed by subsequent processes. 

This scenario has many uncertainties and is hardly 
unique. Major open questions include: in what order 
and on what time scales did scattering of the Excited 
class, accretion of the largest bodies, and depletion of 
the small bodies occur? Were collisions and depletion 
important in the last 3 billion years, with substantial 
cratering, or has the region been largely dormant? Was 
mass loss from the CKB purely due to erosion, or did 
dynamical processes eject significant mass in the form of 
large bodies? What is the present source reservoir for 
Jupiter-family comets? Why is there no "outer" Kuiper 
belt? 

The quantitative data on TNO size and orbital distri- 
butions are improving rapidly. In the next year or two, 
nearly the full available sky will have been surveyed for 



R < 22 TNOs, usin g small telescopes l|Truiillo fc Browiil 
l2003HMoodv et oi.KO O.I.l . Telescopes with effective aper- 
ture « 4 m have surveyed hundreds of square degrees 
ijMillis et aZJl200^ . and will survey thousands of square 
degrees (the CFH Legacy Survey^), or tens of thousands 
(Pan-STARRS®), in the coming 5-10 years, to limits 
of i? < 24, so a full characterization of the accretion- 
dominated size distribution is forthcoming. Exploring 
the erosion-dominated small end is of course more dif- 
ficult. Eight-meter class telescopes can reach i? « 27 
with great effort, and recovery operations are typically 
more expensive than discovery until sky coverage reaches 
tens of square degrees. Wide-area imagers currently com- 
ing into operation on large telescopes will increase the 
number of known i? > 25 TNO orbits past the current 
handful into the tens or hundreds in the next few years. 
Thousands will be possible with a Large Synoptic Sur- 
vey Telescope.^ The situation for further explorations at 
i? > 27 is more bleak: the present study probably repre- 
sents the limits of what can be achieved in this decade, 
unless the HST is tasked to the problem for ^months of 
time, or multi-conjugate adaptive optics becomes capa- 
ble of providing fields of view of many arcminutes. A 
high-throughput space telescope imager such as SNAP^" 
would be capable of improving on the current ACS sam- 
ple by several orders of magnitude circa 2010. The most 
valuable information on D < 50 km TNOs in this decade 
may come from occultati on surveys, e.g., the nearly op- 
erational TAOS project l|Alcock et qiJl200,'^ . And late 
next decade, the New Horizons spacecraft will study the 
composition and structure of the surfaces of Pluto and 
one or more TNOs, providing a cratering and chemical 
record. "'^'^ It is already clear that the joint size-dynamical 
distribution of TNOs is rich in information to constrain 
the physical evolution of the outer Solar System. 
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APPENDIX 
LIKELIHOODS FOR TNO SURVEYS 

Correct analyses of the TNO luminosity function require a proper formulation of the likelihood of a given 

survey detecting a set of objects at magnitudes {m.i}. We cover here some of the subtleties in this process that we have 
glossed over in H3.2I and that have been neglected in other works as well. We then examine how well the simplified 
analyses described in H3 . 21 approximates the likelihoods which might be obtained with full attention to the details, and 
find that, for current sample sizes, these approximations do not affect the inferred bounds on the TNO populations. 

Expected Distribution 

We presume that one is interested in constraining the differential sky density of sources as a function of the 

magnitude R in some standard system. Because the detectability may be a function of some variability parameters 
V, we in fact need to know the joint distribution V) and have 

= J dVJ:{R;V). (Al) 

We assume here that the magnitude R of interest corresponds to the time average over the light curve. The sky density 
will be a function of other variables, e.g. dynamical quantities, but we will leave further generalization to the reader. 

The survey may not be conducted in the standard filter band, so we will denote the mean magnitude in the observed 
filter band as /i, so that we have 



^^{pi-V)= j dRP{^i\R)^{R-V), {A2) 

where we have assumed the object's variability properties to be independent of color. We assume 

Jd^iP{ii\R) = 1. (A3) 

A TNO may be detected by the survey, and after detection it will be assigned an observed magnitude m. The proba- 
bility P(TO|/i; V) of being detected and assigned observed magnitude m is a characteristic of the survey's methodologies, 
which can be determined by proper Monte-Carlo tests. Note that the detection efficiency for a chosen TNO magnitude 
satisfies 

< rjifi-^V) = J dmP{m\n;V) < 1 (A4) 
because the source may not be detected at all. We can likewise integrate over V to define 

P(detection|,) . ^ ^ ^^.3^^' >^ (A5) 
or define a detection probability over R as 

^ Jd,dVr,i,;V)Pi,\R)m-,V)^ 

Note that we can determine rj{ii] V) strictly from analysis of the survey characteristics, because and V fully specify 
the behavior of objects in the imaging. But to estimate 77(/i) or ?7(i?), we also need the distributions of the TNO 
population over variability and/or color, unless the detection efficiency is sufficienctly independent of them. 
The expected distribution of detections from a survey of solid angle is 

g{m)='^^n j dudV P{m\fi;V) J dRP{ii\R)E{R-V). (A7) 

In simple words, the observed distribution g{m) is the convolution of the input distribution S(-R) with the magnitude 
shifts induced by color, variability, measurement error, and detection efficiency. Many real surveys are inhomogeneous, 
with P(m\fi; V) varying with location or time as weather conditions or integration times vary, and in these cases the 
expected g{m) must further be convolved over these variable conditions. 

Likelihoods 

Given the function g(m) derived from some V), we can assign a likelihood to detection of a set {mi, TO2, . . . , ?7ij\/} 
of detections in the survey to be 



L({toJ|E) = dm^ exp 



M 

M ' ' 



dm g{m) 



l[gim,). (A8) 



The form of the likelihood is easily derived by considering each interval 6m of magnitude to be a Poisson process with 
probability SP = g(m) Sm of having one detection and 1 — 6P of zero detections. 
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Comparison to Gladman et al. Form 

IGladman et al\ 1)19981) provide an alternative form of this equation which in our notation and suppressing the 
variabihty parameters, is 

M 

L{{mi}\T,) = exp 



-n J dRr]{R)j:{R) Yi J dR£,{R}J:{R). (A9) 



The fu nction ii{R\m) "describes the uncertainty for the [true] magnitude of object i" according to IGladman et al] 
l)1998j) . but the exact meaning of this is not well understood in the community. It does not mean a simple convolution 
with a Gaussian error term centered on m^. 

The identification of Equation (|A9p with Equation IIASp . up to factors of dm and fl that arc independent of the S 
model, is clear from the following. The exponentiated quantity is in either case N, the expected number of detections, 
given by 

N = j dmg{m) = fl J dm J d^idV P{m\fi;V) J dRP{p\R)Y.{R]V) (AlO) 
= n j dR j dV dmdfi P{m\^;V)P{n\R)^{R]V) (All) 
= nj dRr]{R)i:{R), (A12) 

making use of (|Aip and l)A6|l . 

To equate the terms under the product sign, we note that, if the detection probability is independent of variability, 
we can write 

g{m) j dRP{m\R;V)Y.{R), (A13) 

so that if we identify £i{R) = P{mi\R] V), then Equation (|A9(I becomes equivalent to (|A8(I . Recall that P{m\R; V) is 
the probability of the TNO at R being detected at m, so its integral over m is r]{R) < 1, and there is in any case no 
need to have J dR£i{R) — 1. It is hence not appropriate to take £i{R) to be a Gaussian about with dispersion a 
given by a rough magnitude uncertainty. 

A simplification is possible in the case P{m\R; V) — ri{R)P{m — R), meaning that probability of detection ri{R) is 
independent of the observed magnitude m for given R, and the magnitude measurement error has a fixed distribution 
with unit normalization. This would be the case when the true magnitude is measured in foUowup observations that 
are distinct from the discovery observations. In this case we have 

g{m)= JdRP{m~ R)[nTj{R)]Y.{R) (A14) 
=^ £,{R)=P{m- R)n^s{R). (A15) 

The IGladman et all ((T998|) likelihood can be misinterpreted as assuming that £i{R) = P{mi — R). We see that this 
is incorrect — even in this simplified case it is necessary to have an additional term of rj{R) or iloff(i?) under the R 
integral. In many surveys, however, the detection probability is highly correlated with the observed magnitude m, so 
even this simple form is not applicable. 

Fitting the Amplitude 

As a useful aside, consider the maximum-likelihood estimate of the scale factor s that normalizes a candidate TNO 
distribution function of the form 

Y.{R:V) = sY.o{R\V). (A16) 

The expected distribution will clearly scale as g{m) = sgo{m), with the obvious notation. The likelihood (|A8|I then 
has all of its dependence on s in terms 

Locs^^e-^^", (A17) 

where No is the expected number of detections N from Equation (|A10(I for E = So. If we adopt a logarithmic prior 
p{s) oc s~^, then the Bayesian expectation values and probability distribution for s are 

s={s)^M/No (A18) 

(N) = M (A19) 

pis> So)=PiM,SoNo), (A20) 
where the last right-hand side is the incomplete gamma function. 

For M = (no detections), the logarithmic prior leads to unnormalizable posterior distributions, and we revert to a uniform prior. 
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If we are interested only in the shape Eq of the TNO distribution, we will marginalize over the nuisance parameter 
s, which for the logarithmic prior yields 



L({to,}|So)- j dsp{s)L{{m^}\^o,s) (A21) 
= 7V-*^r(Af)[].9o(m,) (A22) 

i 

= M-^'^r(Af) Y[ s go{m,) (A23) 

i 

= Af-*^r(M)e*^L({m,}|Eo,s). (A24) 
This analytic marginalization speeds up our likelihood analyses for the shape of S. Note that the formulae fail for 
M = 0, which is not surprising, since we cannot constrain the shape of S if we have no detections. 

Detection Cutoff 

In the following we derive some exact results for simple cases, and then examine the error made in simple approxi- 
mations. We will assume that in the neighborhood of the survey limit, the source density behaves as = SqIO"''^. 

The simplest case to analyze is when the survey assigns a detection magnitude pi to each TNO, and those brighter 
than a cutoff po are accepted into the survey. From measurement noise alone, the detection magnitude usually has 
a Gaussian distribution about the true magnitude R with a standard deviation of cf{R) mag. If color and variability 
are important, we can consider these as additional sources of scatter in the P{p\R) distribution, as quantified in later 
sections, for now we consider measurement noise alone. 

Consider two cases: in the first, which we'll call a snapshot survey, the survey is halted after its detection observations, 
so that riii = Pi. If we further assume that <j{R) varies little near the completeness limit, then we have 

g{m) = f^SolO"™ exp[a2 ln2(10) Q{po ~ m) (A25) 
ov 

N=—^WP'-exp[a^\n\lO)a^/2] (A26) 

r){R) = (27r(T2)-i/2 P exp[-(i? - p)y2a^] (A27) 

^ — OO 

= ^erk[{R- Po)/V2a]. (A28) 

Here O is the step function. In this case, ignoring the diff erence between observed and detected magnitudes — 
that is, using instead of g{m) in the likelihood IjASp for the {m^} — amounts to a normalization error of 

exp[Q;^ln (10)cr^/2]. The completeness follows Equation Q with w = ^/2a. 

In a second type of survey, which we'll call a magnitude followup survey, the detected objects are followed up with 
more integration, so that the errors on magnitude are substantially reduced, and we end up with = Ri, in the 
absence of color terms. Our ACS survey is well approximated in this way, since the faintest detected object has 
S/N > 20 in the full dataset, hence a magnitude uncertainty (in F606W) of < 0.05 mag. In this case we have 

g{m) = nj:{RMR), (A29) 
where the detection efficiency ri{R) is from Equation IIA28I) . and is in fact independent of the form of We note 

that in the case of a magnitude- followup survey, the approximate likelihood in Equation Q is exactly correct. 

Size of the Errors 

Suppose we were to conduct a snapshot survey that detects objects at { mj}, and also a Monte Carlo survey to 
determine the form of ri{R), but then we get lazy in implementing Equation (jA8|l and ignore the difference between 
m and R in calculating our likelihoods, using 

g{m) = nY.{m)Q{po - m) (A30) 
instead of the proper form in Equation (|A25ll . Then we see that we are in effect using a luminosity function that is in 
error by the normalization factor F — exp[Q!^ In^(lO) cr^/2]. We can bound this factor by noting that a TNO detection 
must have S/N > 6 to be reliably distinguished from noise in most detection schemes, so we expect the noise on the 
detected magnitude to be cr < 2.5/[ln(10) (S/N)] < 0.18 mag. In this case we have F 1 + 0.04(a/0.7)^. First we 
note that this bias factor, if constant across all survey data, only affects the normalization of S and does not affect 
our conclusions about the shape of S. For \a\ < 0.7, the resultant 4% bias would be smaller than the Poisson errors in 
our analysis, which contains only 130 total objects and is hence incapable of detecting < 10% shifts in the presumed 
at any magnitude. 

If the bias factor F varies with R than we could infer an incorrect shape for E, but again with our sample of 
130 TNOs the difference is below the Poisson noise. Note that we conclude \a\ < 0.7 for all samples and magnitudes, 
save perhaps the bright end of the CKBO population, where we could incur a bias somewhat above 10% — but the 
Poisson noise is also higher in this subsample, so again the bias is negligible. Note that for surveys with significant 
magnitude followup, there are no approximations in taking g{m) = r2ry(m)E(m), so there is no bias at all. This is 
the case for our ACS data, and also for the La and TB data, since we have accurate followup mean R magnitudes for 
nearly all detections. 
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Inhomogeneous Surveys 

If we take this simple case literally, we would expect all surveys to report a completeness function ri{R) that follows 
Equation (|2Jl with w « {S/N)~^ /^/2. This is true in our case, with a limiting S/N = 8.2 and w = 0.08. But most 
surveys report completeness functions ri{R) that are significantly broader than this form. This is generally due to 
inhomogeneity, as ground-based surveys must deal with varying seeing and coverage, etc., that the ACS survey does 
not. In effect the survey area f2 is divided into sub- areas U,k that have cutoffs pk- The expected distributions gk{m) 
for each sub-survey may be summed to yield the total g{m) for the survey. It is convenient to define the total survey 
area which reaches a chosen threshold: 

np{p) = ^fc- (^31) 
Pk>p 

For a snapshot survey, we obtain: 

5(m) = rap(m)SolO"™ (A32) 

^(i?) = (27ra2)- 1/2 jdp^^ cxp[-(i? - pY/2a^] {A33) 

= JaJ ^^^7 "'^"[^^ ~ P^^^''^- ("^^^^ 

The falloff in completeness ri{R) is now the convolution of the inhomegeneity with the measurement error, so will be 
more extended than either. But the bias F is not affected by the inhomogeneity, only by the measurement error, so 
we see that the approximation in Equation ljA30|l can be replaced with 

g{m) « rjp(TO)E(m) (A35) 

with equal accuracy, sufficient for our purposes. 
For an inhomogeneous magnitude-followup survey, we find that once again Equation ((HJ is exact. 

Color Errors 

The comparison of the observed-band magnitude p with a density function E defined over a standard magnitude 
R can be thought of as a convolution of the expected S(i?) with the color distribution P{p\R). We have assumed a 
single value for our i^606W^ — R color conversion, but there is a spread in colors which may be viewed as an additional 
source of stochastic uncertainty on the R magnitude of our detections. The sta ndard deviation oi V — R colors is 
O.IO mag (jTegler fc Romani shin"2003 HTegler. Romanishin. fc C onsolm agndl2003(l . Since F&O&W is nearly the union 
of V and i?, similar to the "VT?" filter used by several ground-based surveys, we estimate the scatter in FQOQW — R 
to be avR ~ 0.05 mag. 

This scatter may be treated as a measurement error on R that is added in quadrature to the measurement noise. 
The < 0.1 mag uncertainty due to the color conversion should have negligible effect on the interpretation of at 
least for the surveys under current consideration which have < 100 detections and are in well-calibrated filters bounded 
by the V and R passbands. Note, however, that the two brightest surveys, TB and La, have very poorly defined filter 
bands, so the color uncertainties could be larger here. For most of these bright objects, however, a high-precision R 
magnitude has been found in the literature, so the color or photometry uncertainties are minimal. 

Variability 

The variability of the TNOs can influence both the detectability and the observed magnitude nrii that is assigned to 
a discovery. The effect of variability on the detection function P{m\p\ V) has not been accurately quantified for any 
published survey, and furthermore the variability dependence of S(-R, V) is very poorly constrained, so at this point 
it is basically impossible to meaningfully incorporate variability into the likelihood functions. 

We can nonetheless show with some basic modelling that the effect of variability on current analyses should be 
minimal given the present Poission uncertainties. The true form of P{m\p] V) for real-life surveys is probably rather 
complex: most surveys require detection in two or more observation sets seperated by hours to days, and the known 
variable TNOs have periods of 3-13 hours. We will consider some simplified forms of variability selection. 

Snapshot Survey 

A "snapshot" survey in this context observes the TNO only for a time short compared to the lightcurve period, so 
the observed magnitude mi is an instantaneous sample of the light curve, and the detection probability is a function 
solely of rui. We ignore color effects and measurement noise here, for simplicity. For the snapshot survey, let R{t) be 
the magnitude at time t. We then have 

P{m\R; V) = e(m - po)P{m\R{t))P{R{t)\R] V), (A36) 

where is a step function. P{R{t)\R; V) is simply the distribution of instantaneous magnitudes for the chosen light 
curve. So g{m) is the convolution of the intrinsic distribution S(-R) with the distributions due to variability and 
measurement noise. The variability in this case just acts to broaden the detection cutoff vs mean magnitude, just like 
a source of measurement noise. 
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How broad is the variability kernel? For a sinusoidal light curve with peak-to-peak amplitude A, the distribution 
P(m - R) oc [1~ {m~ Rf j^A^Y^I'^. When a power -law is convolved with this distribution, the resultant 

= AN I dm is enhanced by the factor 

F = 1 + (Aa In 10) Vl6 + (Aa In 10)'*/1024 + ... (A37) 

Even in the worst case of a square-wave light curve, the bias factor is F = cosh(AQ; ln(10)/2), which is « 2 times 
larger. 

Our 2003 BFgi, along with 2 001 QG2fl8, are the most highly variable of t he ~ 60 TNOs tested for variability to 
date, with A « 1.1 mag (^TriUing fc Bernsteinll200l iSheopard fc Jewitlil2004j) . With this A and a < 0.7, the surface 
density for instantaneous magnitudes is F < 1.20 higher than the S(-R) for mean magnitudes. 

One might expect small (faint) TNOs to be less symmetric and hence more variable than bright TNOs, though this is 
not yet borne out. But it can never be the case that they have uniformly high amplitudes A, because pole-on rotators 
must have zero amplitude. The formulae of ' Lacerda fc Luu ( 2003) for the light-curve amplitudes of triaxial ellipsoids 
in minimum-energy rotation can be used to show that a population of identically-shaped TNOs with randomly oriented 
rotation axes will have a roughly uniform distribution of light-curve amplitudes between and Aj^ax- We might expect 
a similar result for variability due to surface features rather than gross body shape. In this case the population-averaged 
enhancement factor will be 

F « 1 + (Anaxo In 10) V48. (A38) 

Take, as a pessimistic case, a population of TNOs that arc all sufficiently asymmetric to have light-curve amplitudes 
as large as 2003 BFgi if viewed 90° from the rotation pole. Then we find the difference between the instantaneous- 
and mean-magnitude distributions to be F < 1.07, less than 7%. Thus the error made by failing to distinguish these 
cases is below the Poisson noise of the current samples. 

Worst- Case Models 

The worst possible bias in the TNO sample due to variability would be if the selection were based on peak (or 
minimum) flux over the light curve. A simplistic scenario that might spawn such a selection effect is for a survey that 
consists of many repeated observations of each target, and a "detection" consists of the TNO rising above the threshold 
in one or more of the epochs. Once detected, the TNO is found in all epochs, and we report a mean magnitude. In this 
scenario, the reported magnitude is an unbiased estimator of the true R, but the detection criterion is a step function 
of the peak flux over the light curve, ri{R, V) = 0[pc — {R — A/2)], for a light curve with peak-to-peak amplitude A. 
Note that in this case, as for the snapshot model in the previous paragraph, there is a bias toward detection of highly 
variable TNOs. Hence increased variability in the faint population would tend to steepen the derived S(i?), as opposed 
to the flattening we detect in the ACS sample. 

Conversely imagine that the detection criterion is that the TNO be above threshold in all of the epochs. In this case 
r]{R, V) will be a function solely of the minimum flux over the light curve, and the detection is biased against high 
variability at a given mean R. 

Again take a pessimistic view that A is uniformly distributed between and Amax ~ 1 mag. Then (in the absence 
of noise) the selection function r]{R) drops linearly from unity at po — A/2 to zero at po- This closely resembles the 
error-function form in Equation Q with mso = po — ^max/4 and w = A,„ax/8. Hence in this worst-case scenario, the 
effect of variability is to shift our completeness limit brightward by < 0.3 mag and broaden the cutoff by < 0.15 mag. 

Is a Simple Treatment Sufficient? 

In analysing the surveys taken to date, we have taken a simplified form for the likelihood, embodied by Equation 
We have seen that this form is exact in the case of a magnitude-foUowup survey with a fully realistic Monte-Carlo 
derivation of ficif (i?) = ftp{R) that is properly marginalized over the variability and color of the TNOs. Technically 
the simplified likelihood is incorrect for our analysis because: 

• The middle-magnitude surveys have no foUowup and are "snapshot" surveys, meaning that g{m) from Equa- 
tion ljA32|l is more appropriate. We have shown that the F factor is smaller than Poisson er rors, a nd the difference 
between f2p(m) and r]{m) is only a convolution by the relatively sharp function Equation (|A28(I from magnitude 
errors of cr < 0.2 mag. If we avoid the regions where rj < 0.2, then the effect on the likelihood is small. 

• The Monte-Carlo estimates of r]{R) in most surveys ignore errors in color terms. The effect is to broaden ri{R) 
by < 0.05 mag in quadrature. 

• The Monte-Carlo estimates of completion ignore variability. For short-exposure surveys this is equivalent to a 
< 7% mis-measure of the underlying I](i?). In the worst imaginable case of selection effects, it could lead to a 
shift in the completeness limit by < 0.3 mag and a broadening of < 0.15 mag. 

• For technical reasons the Monte-Carlo for the ACS data yields ftp{rn) rather than the slightly broader ri{R) 
function. 

• Some readers may not trust the completeness tests or magnitude error estimates of our or other surveys. 
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To test the effect of these approximations on our resuhs, we have re-analysed the survey data after artificiaUy broad- 
ening the »7(i?) cutoff width w by 0.3 mag for the ACS and/or TB surveys. We also shift the cutoff mso of the ACS 
completion function brightward by 0.3 mag to test the import of a substantial mis-estimate of our completeness or 
variability treatments. 

Note that we do not apply these completeness changes to all the surveys, just the ones at the ends of the magnitude 
ranges, because we are more interested in errors that might have influenced the shape of the fits. We find that 

no combination of these tweaks to the assumed r]{R) functions causes a significant shift in the {011,02} contours of 
FigureElor in any of the quantities calculated herein. When we shift m^o of the ACS survey brightward, the ability to 
discriminate 02 for the two classes is weakened (distinct 02 values are still strongly suggested) , and the extrapolation 
to the JFC source population can be higher. Our main conclusions, however, remain valid: a break in the power- 
law S(m) is required, there is a significant deficit of small bodies, the CKB and Excited classes have distinct size 
distributions, with the latter having clearly shallower ai. We therefore conclude that a simplified treatment of the 
likelihoods is currently adequate. Simply put, the ACS survey has only 3 detections, so there is no need to be too 
careful with subtle variability or color effects. 

When uniform TNO surveys with 3> 100 detections are available, however, a more careful treatment will be required. 
In particular, careful attention to the effect of variability may be needed, and the likelihood analyses should perform 
the proper integrals to get g{m) in the case that accurate foUowup magnitudes are not available, because the distinction 
between the discovered magnitude m and the true magnitude R is important at the few-percent level. 



